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Abstract 

We study the spectral stability of a family of periodic wave trains of the Korteweg-de 
Vries/Kuramoto-Sivashinsky equation dtv + vdxV + d^v + S (d^v + d^v) ~ 0, S > 0, in 
the Korteweg-de Vries limit (5 — >■ 0, a canonical limit describing small-amplitude weakly 
unstable thin film flow. More precisely, we carry out a rigorous singular perturbation 
analysis reducing the problem to the evaluation for each Bloch parameter ^ S [0, 27r] of 
certain elliptic integrals derived formally (on an incomplete set of frequencies/Bloch pa- 
rameters, hence as necessary conditions for stability) and numerically evaluated by Bar 
and Nepomnyashchy [BN], thus obtaining, up to machine error, complete conclusions 
about stability. The main technical difficulty is in treating the large-frequency and 
small Bloch-parameter regimes not studied by Bar and Nepomnyashchy [BN], which 
requires techniques rather different from classical Fenichel-type analysis. The passage 
from small-(5 to small-^ behavior is particularly interesting, using in an essential way an 
analogy with hyperbolic relaxation at the level of the Whitham modulation equations. 
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1 Introduction 

In this paper, we study the spectral stability of periodic wave trains of the Korteweg-de 
Vries/Kuramoto-Sivashinsky (KdV-KS) equation 

(1.1) dtu + ud^u + dlu + S{dlu + d^u) =0, Vi>0,VxGM, 

with < S <^ 1. When 6 = 0, equation (1.1) reduces to the well-studied Korteweg-dc Vries 
(KdV) equation, which is an example of a completely integrable infinite dimensional Hamil- 
tonian system. As such, the KdV equation is solvable by the inverse scattering transform, 
and serves as a canonical integrable equation in mathematical physics and applied mathe- 
matics describing weakly nonlinear dynamics of long one-dimensional waves propagating in 
a dispersive medium. 

When (5 > on the other hand, equation (1.1) accounts for both dissipation and disper- 
sion in the medium. In particular, for < (5 ^ 1 it is known to model a thin layer of viscous 
fluid flowing down an incline, in which case it can be derived either from the shallow water 
equations 

dth + dx{hu) = 0, dt{hu) + dx{hu^ + ^^) = h-u^ + udxikd^u). 

as -F — 7- 2"'" {F being the Froudc number, with F = 2 the critical value above which steady 
constant-height flows are unstable) or from the full Navier-Stokes equations if < i? — i?c <C 
1 {R c being the critical Reynolds number above which steady Nusselt flows are unstable) in 
the small amplitude/large scale regime: see [Wi, YY] for more details. For other values of 



1 INTRODUCTION 



3 



S, (1.1) serves as a canonical model for pattern formation that has been used to describe, 
variously, plasma instabilities, flame front propagation, or turbulence/transition to chaos in 
reaction-diffusion systems [SI, S2, SM, K, KT]. 

Here, our goal is to analyze the spectral stability of periodic traveling wave solutions 
of (1.1) with respect to small localized perturbations in the singular limit ^ ^ 0. In this 
limit the governing equation (1.1) may be regarded as a dissipative (singular) perturbation 
of the KdV equation, for which it is known that all periodic traveling waves are spectrally 
stable to small localized perturbations; see [BD, KSF, Sp]. However, as the limiting KdV 
equation is time-reversible (Hamiltonian) , this stability is of "neutral" (neither growing nor 
decaying) type, and so it is not immediately clear whether the stability of these limiting 
waves carries over to stability of "nearby" waves in the flow induced by (1.1) for \6\ <^ 1. 
Indeed, we shall see that, for different parameters, neutrally stable periodic KdV waves may 
perturb to either stable or unstable periodic KdV-KS waves, depending on the results of a 
rather delicate perturbation analysis. 

Our analysis, mathematically speaking, falls in the context of perturbed integrable sys- 
tems, a topic of independent interest. In this regard, it seems worthwhile to mention that 
the proof of stability of limiting KdV waves, and in particular the explicit determination 
of eigenvalues and eigenfunctions of associated linearized operators on which the present 
analysis is based, is itself a substantial problem that remained for a long time unsettled. 
Indeed, by an odd oincidence, both the original proof of spectral stability in [KSF, Sp] and 
a more recent proof of spectral and linearized stability in [BD] (see also the restricted non- 
linear stability result [DK]) were accompanied by claims appearing at about the same time 
of instability of these waves, an example of history repeating itself and indirect indication 
of the difficulty of this problem. 

However, our motivations for studying this problem come very much from the physical 
applications to thin-film flow, and particularly the interesting mctastability phenomena 
described in [PSU, BJRZ, BJNRZ3] (see Section 1.1 below). Interestingly, our resolution of 
the most difficult aspect of this problem, the analysis of the small- Floquet number/small-5 
regime, is likewise motivated by the associated physics, in particular, by the formal Whitham 
equations expected to govern long-wave perturbations of KdV waves, and an extended 
relaxation- type system formally governing the associated small-5 (KdV-KS) problem. 

The identification of this structure, and the merging of integrable system techniques with 
asymptotic ODE techniques introduced recently in, e.g., [JZ2, PZ, HLZ, BHZ] (specifically, 
in our analysis of frequencies |A| G [C, C(5~^]), we regard as interesting contributions to 
the general theory that may be of use in related problems involving perturbed integrable 
systems. Our main contribution, though, is to the theory of thin-film flow, for which the 
singular limit 6—^0 appears to be the canonical problem directing asymptotic behavior. 

We begin by defining the notion of spectral stability of periodic waves of the 4th-order 
parabolic system (KdV-KS), following [BJNRZl], as satisfaction of the following collection 
of nondcgcneracy and spectral assumptions: 

• (HI) The map ii" : ^ taking {X,b,c,q) ^ {u,u' ,u"){X,b,c,q) - bis full rank 
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at {X,b,c,q) G H ^({0}) where {u,u' ,u"){-;b,c,q) is the solution of 

5(u"' + u) + u" + (u — c)u = q, {u, u ,u"){0; b, c, q) = b. 

• (Dl) aL^R){L) C {A G C I SR(A) < 0} U {0}, where 

L = -6 {d^^ + dl)-dl-dAu-c) 
denotes the hnearized operator obtained by Hnearizing (1.1) about u = u{.;b,c,q). 

• (D2) o-i2^^(o,x)(L^) C {A G C|SRA < for some ^ > and any ^ G [-7r/X,7r/X) 
where 

Ldf] = -S {(dx + iO^f + (5. + - {d. + - {d. + ^6 {{u - c)f) 

denotes the associated Bloch operator with Bloch-frequency ^. 

• (D3) A = is an eigenvalue of the Bloch operatorLo of algebraic multiplicity two. 

Under assumptions (HI), {D3), standard spectral perturbation theory implies the exis- 
tence of two eigenvalues Aj(^) G o'(L^) bifurcating from the (A. ^) = (0,0) state of the form 
"^jiO = ^'^iC + o(C)- Assumption (Dl) ensures that aj G M. To ensure analyticity in ^ of 
the critical curves A^-, we assume further: 

• (H2) The coefficients aj are distinct. 

The above definition of spectral stability is justified by the results of [BJNRZl], which 
state that, under assumptions (H1),(H2) and (D1)-(D3), the underlying wave u is L-'^(M) fl 
H*{R) L°°(M) nonlinearly stable; moreover, if u is any other solution of (1.1) with data 
sufficiently close to u in L-'^(M)ni!f^(M), for some appropriately prescribed i/j„ the modulated 
solution u{- — ip{-,t),t) converges to u in L^'(M), p G [2,oo]. 

This is to be contrasted with the notion of spectral stability of periodic waves of Hamilto- 
nian systems, which, up to genericity conditions analogous to (H1)-(H2) and (D3), amounts 
to the condition that the associated linearized operator analogous to L have purely imagi- 
nary spectrum. That is, in order that a (neutrally) stable periodic wave of (KdV) perturb 
under small 5 > to a stable periodic wave of (KdV-KS), its spectra must perturb from 
the imaginary axis into the stable (negative real part) complex half-plane 

The main goal of this paper, therefore is to establish by rigorous singular perturbation 
theory a simple numerical condition guaranteeing the existence of periodic traveling wave 
solutions of (1.1) satisfying the above hypotheses (H1)-(H2) and (D1)-(D3) for sufficiently 
small (5 > 0: more precisely, determining whether the neutrally stable periodic solutions of 
(KdV) perturb for small 5 > to stable or to unstable solutions of (KdV-KS). 

Remark 1.1. The methods used in [BJNRZl] to treat the dissipative case 6 > 0, based 
on linearized decay estimates and variation of constants, are quite different from those 
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typically used to show stability in the Hamiltonian case. The latter are typically based on 
Arnolds method, which consists of finding sufficiently many additional constants of motion, 
or "Casimirs, " that the relative Hamiltonian becomes positive or negative definite subject to 
these constraints (hence controlling the norm of perturbations); that is, additional constants 
of motion are used to effectively "excise" unstable (stable) eigenmodes of the second variation 
of the Hamiltonian. This approach is used in [DK] to show stability with respect to nX- 
periodic perturbations for arbitrary G N, where X denotes the period of the underlying 
KdV wave train. However, as spectra in the periodic case is purely essential, such 
unstable (stable) eigenmodes are uncountably many, and so it is unclear how to carry out 
this approach for general perturbations. Indeed, to our knowledge, the problem of stability 
of periodic KdV waves with respect to general H'^ perturbations remains open. 

The spectral stability of periodic wave-train solutions of (1.1) itself has a long and 
interesting history of numerical and formal investigations. In [CDK], the authors studied 
numerically the spectral stability of periodic wave trains of 

dtu + udxU + jd^u + d'^u + d^u = 0, 

which is, up to a rescaling, equation (1.1) with S = and showed the stabilizing effect 
of strong dispersion (large 7/small ds). As 7 is increased from to 8 ((5 G (1/8, 00)), only 
one family of periodic waves of Kuramoto-Sivashinsky equation survives and its domain of 
stability becomes larger and larger and seems to "converge" to a finite range (Li,L2) with 
Li = 27r/0.74 and L2 ~ 27r/0.24. In [BN], Bar and Nepomnyashchy studied formally the 
spectral stability of periodic wave trains of (1.1) as (5 — )■ 0. More precisely, for a fixed Bloch 
wave number ^ the associated eigenvalues X{C^S) G cr(L^) are formally expanded as 

(1-2) X{5,0 = XoiO + S\iiO + OiS^), 

where Ao(^) G Ri is an eigenvalue associated to the stability of periodic waves of KdV, 
known explicitly (see [BD]), and Ai(^) is described in terms of elliptic integrals. Then, the 
authors verified numerically, using high-precision computations in MATHEMATICA (see 
Appendix B, [BN]), that inax^g[_^/x ,r/x) ^('^i(O) < 0' consistent with stability, on the 
band of periods X G (Li, L2) with Li ~ 8.49 and L2 ~ 26.17, which arc approximately the 
bounds found in [CDK]. Similar bounds were found numerically in [BJNRZl] by completely 
different, direct Evans function, methods, with excellent agreement to those of [BN]. 

However, the study of Bar and Nepomnyashchy is only formal and in particular, as 
mentioned in [BN], it is not valid in the neighbourhood of the origin (A,^) = (0,0). In 
particular, as we show in Section 4, the description (1.2) is valid only for < ||t <C 1 hence, 
for any given 5 > it is not possible to conclude from this expansion spectral stability of 
an associated periodic wave train of (1.1). (The formal derivation of this expansion in [BN] 
is valid for |(A,^)|cxM bounded from zero.) 

Likewise, the numerical study in Section 2 of [BJNRZl], which is not a singular per- 
turbation analysis, but rather a high-precision computation down to small but positive 5,^ 

^Minimum value 5 = .05, as compared to <5 = .125 in [CDK]; see Table 3, [BJNRZl]. 
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gives information about 5^0 only at finite scales, hence in effect omits an 0{6 x TOL) 
neighborhood of the origin, where TOL is the relative precision of the computation. Thus, 
though very suggestive, neither of these computations gives conclusive results about stabil- 
ity in the (5 — >■ limit, and, in particular, behavior on a 5 -neighborhood of the origin is not 
(either formally or numerically) described. 

In this paper, we both make rigorous the formal singular perturbation analysis that was 
done in [BN] and extend it to the frequency regimes that were omitted in [BN] , completing 
the study of the spectrum at the origin and in the high-frequency regime. More precisely, 
we carry out a rigorous singular perturbation analysis reducing the problem to the study 
of Bloch parameters ^ G [C(5, 27r] (see Section 3 for definition of Bloch parameter) and 
eigenvalues |3?A| < C(5, |S^A| < C, C > sufficiently large, on which the computations of 
[BN] may be justified by standard Fenichel-type theory. 

The exclusion of high frequencies is accomplished by a standard parabolic energy esti- 
mate restricting |3A| < C5^^/'^ followed by a second energy estimate on a reduced "slow," 
or "KdV," block restricting |9A| < C; see Lemma 3.1 and Proposition 3.3 in Section 3.1 
below. For related singular perturbation analyses using this technique of successive reduc- 
tion and estimation, see for example [MZ, PZ, Z, JZ2] and especially [BHZ], Section 4. 
The treatment of small frequencies proceeds as usual by quite different techniques involv- 
ing rather the isolation of "slow modes" connected with formal modulation and large-time 
asymptotic behavior. 

At a technical level, this latter task appears quite daunting, being a two-parameter 

bifurcation problem emanating from a triple root A = of the Evans function at ^ = 
(5 = 0, where the Evans function -^(A, ^, 5) (defined in (3.4) below) is an analytic function 
whose roots A for fixed (5 > and ^ G [—tt/X,t:/X) comprise the L^-spectrum of the 
linearized operator about the periodic solution. However, using the special structure of the 
problem, we are able to avoid the analysis of presumably complicated behavior on the main 
"transition regime" A := {((5,0 I ^ l?l/<^ ^ C*}, C 2> 1 and only examine the two 
limits \^\I5 — > 0, +oo on which the problem reduces to a pair of manageable one-parameter 
bifurcation problems of familiar types. 

Specifically, wc show that (small) roots A of the Evans function cannot cross the imag- 
inary axis within A, so that stability properties need only be assessed on the closure of 
A-complement, with the results then propagating by continuity from the boundary of A 
to its interior. This has the further implication that stability properties on the wedges 
<^ 1 and \(,\/5 ^ 1 are linked (through A), and so it suffices to check stability on 
the single wedge \S,\/5 ^ 1, where the analysis reduces to computations carried out in 
[BN]. Indeed, the situation is simpler still: stability on the entire region ^ 1 reduces 
by the above considerations to validity of a certain "subcharacteristic condition" relating 
characteristics of the Whitham modulation equations for KdV (the limit |^|/<5 — >■ oo) and 
characteristics of a limiting reduced system as () — > (the limit 0). 

As the above terminology suggests, there is a strong analogy in the regime |^|,(5 <C 1 to 
the situation of symmetric hyperbolic relaxation systems and stability of constant solutions 
in the large time or small relaxation parameter limit [SK, Yo, Ze], for which a similar 
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"noncrossing" principle reduces the question of stability to checking of Kawashima's genuine 
coupling condition, which in simple cases reduces to the subcharacteristic condition that 
characteristics of relaxation and relaxed system interlace. Indeed, at the level of Whitham 
modulation equations, the limit as \C\/S oo may be expressed as a relaxation from the 
Whitham modulation equations for KdV to the Whitham equations for fixed 6 in the limit as 
(5 —7- 0, a relation which illuminates both the role/meaning of the subcharacteristic condition 
and the relation between KdV and perturbed systems at the level of asymptotic behavior. 
These issues, which we regard as some of the most interesting and important observations 
of the paper, are discussed in Section 4.2. 

The final outcome, and the main result of this paper, is that stability- whether spectral, 
linear, or nonlinear- of periodic solutions of (1.1) in the KdV limit (5—7-0 apart from the 
periods Xi Ri 8.49 and X2 ~ 26.17 at which Ind(X) = max_jr/x<5<7r/x 3ftAi(^) precisely 
vanishes, is determined by d-{X) := sgnlnd(X), or, equivalently, the values of the elliptic 
integrals derived by Bar and Nepomnyashchy in [BN], with <t(X) < corresponding to 
stability. 

1.1 Discussion and open problems 

As we have emphasized above, our analysis, though convincing we feel, does not constitute a 
numerical proof, but rather a "numerical demonstration," in the sense that the computations 
of [BN] on which we ultimately rely for evaluation of d'{X) are carried out with high precision 
and great numerical care, but not with interval arithmetic in a manner yielding guaranteed 
accuracy. However, there is no reason that such an analysis could not be carried out- we 
point for example to the computations of [M] in the related context of stability of radial 
KdV-KS waves- and, given the fundamental nature of the problem, this seems an important 
open problem for further investigation. 

Indeed, numerical proof of stability or instability for arbitrary nonzero values of 6, 
verifying the numerical conclusions of [BJNRZl], or of Evans computations in general, 
though considerably more involved, seems also feasible, and another important direction 
for future investigations. 

The particular limit S ^ studied here has special importance, we find, as a canonical 
limit that serves (as discussed at the beginning of the introduction) as an organizing center 
for other situations/types of models as well, and it has indeed been much studied; sec, for 
example, [EMR, BN, PSU], and references therein. As discussed in [PSU, BJNRZ3, BJRZ], 
it is also prototypical of the interesting and somewhat surprising behavior of inclined thin 
film flows that solutions often organize time-asymptotically into arrays of "near-solitary 
wave" pulses, despite the fact that individual solitary waves, since their endstates necessarily 
induce unstable essential spectrum,^ are clearly unstable. 

To pursue the analogy between modulational behavior and solutions of hyperbolic- 
parabolic conservation or balance laws that has emerged in [OZ, Se, BJNRZl, BJNRZ2], 
etc., and, indeed, through the earlier studies of [FST] or the still earlier work of Whitham 
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[W], we feel that the KdV hmit S of (1.1) plays a role for small-amplitude periodic 
inclined thin film flow analogous to that played by Burgers equation for small-amplitude 
shock waves of general systems of hyperbolic-parabolic conservation laws, and the current 
analysis a role analogous to that of Goodman's analysis in [Gol, Go2] of spectral stability 
of general small-amplitude shock waves by singular perturbation of Burgers shocks.^ 

The difference from the shock wave case is that, whereas, up to Galilean and scaling in- 
variances, the Burgers shock profile is unique, there exists up to invariances a one-parameter 
family of periodic waves of KdV, indexed by the period X, of which only a certain range are 
stable. Moreover, whereas the Burgers shock profile is described by a simple tanh function, 
periodic KdV waves are described by a more involved parametrization involving elliptic 
functions. Thus, the study of periodic waves is inherently more complicated, simply by 
virtue of the number of cases that must be considered, and the complexity of the waves in- 
volved. Indeed, in contrast to the essentially geometric proof of Goodman for shock waves, 
we here find it necessary to use in essential ways certain exact computations coming from 
the integr ability /inverse scattering formalism of the underlying KdV equation. 

Plan of the paper. In Section 2, we compute an expansion of periodic waves of KdV- 
KS in the limit 5 — )■ by using Fenichel singular perturbation theory. In Section 3, we 
justify rigorously the formal spectral analysis in [BN]: we provide, first, a priori estimates 
on the size of unstable eigenvalue and show that they arc necessarily of order 0(1) as 5 — t- 0. 
Then we compute an expansion of both the Evans function and eigenvalues with respect to 
6. This analysis holds true except in a neighborhood of the origin from which spectral curves 
bifurcate. In Section 4, we compute the spectral curves in the neighborhood of the origin 
and show that spectral stability is related to subcharacteristic conditions for a Whitham's 
modulation system of relaxation type. 

2 Expansion of periodic traveling-waves in the KdV limit 

For < (5 <C 1, equation (1.1) is a singular perturbation of the Korteweg-de Vries equation 
(2.1) dtu + ud^u + dlu = 0, 

where the periodic traveling wave solutions may be described with the help of the Jacobi 
elliptic functions. In [EMR], periodic traveling wave solutions of (1.1) are proved to be 
5-close to periodic traveling wave solutions of (2.1) and, furthermore, an expansion of 
these solutions with respect to 5 is found. We begin our analysis by briefly recalling the 
details of this expansion. Notice that (1.1) admits traveling wave solutions of the form 
u{x,t) = U{x — ct) provided the profile U satisfies the equation 

{U - c)U' + U'" + 5{U" + U"") = 0, 



^See also the related [PZ, PreS], more in the spirit of the present analysis. 
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where here ' denotes differentiation with respect to the traveling variable oj = x — ct. Due 
to the conservative nature of (1.1) this profile equation may be integrated once yielding 

(2.2) ^-cU + U" + S{U' + U"') = q, 

where g G M is a constant of integration. By introducing x = U,y = U' and z = U" + U, 
we may write (2.2) as the equivalent first order system 

(2.3) x' = y, y' = z-x, Sz' = —z + q + {c + l)x — —. 
Setting (5 = in (2.3) yields the slow system 

x^ 

(2.4) 2; = g + (c+ l)x - y, x' = y, y' = q + cx-—, 

which is equivalent to the planar, integrablc system governing the traveling wave profiles for 
the KdV equation (2.1). Utilizing the well-known Fcnichel theorems, we are able to justify 
the reduction and continue the resulting KdV profiles for < (5 ^ 1. To this end, we define 



Mo = <^ {x,y,z) e 



x^ 

z = q+{c+l)x - — =: foix) 

and recognize this as the slow manifold associated to (2.3). It is readily checked that this 
manifold is normally hyperbolic attractive, and so a standard application of the Fenichel 
theorems yields the following proposition. 

Proposition 2.1. For 5 > sufficiently small, there exists a slow manifold Ms invariant 
under the flow of (2.3) that is written as 

M5 = {(x,2/,z) gR3 I z = fs{x,y)], fs{x,y) = h{x) + 5fi{x,y) + 5^ f2{x,y) + 0{5^), 



J.2 

fi{x,y) =xy- (c+ l)y, f2{x,y) = -y^ + {c + l- x){q + cx- y). 

The expansion of fs is obtained by inserting z = fsix,y) into (2.3) and identifying the 
powers in 6. Then by plugging this expansion into (2.3)2, one finds the reduced planar 
system: 

(2.5) x' = y, y' = q + cx- — + S{xy-{c+l)y) + 0{S'^), 
or equivalently the scalar equation 

(2.6) x" = q + cx - — + 6{x-c- l)x' + 0{S'^). 

Now, we seek an asymptotic expansion of the solutions of (2.5) in the limit (5—7-0. An 
easy way of doing these computations to any order with respect to S is to follow the formal 
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computations in [BN], which are now justified here with Fenichel's theorems. To begin, 
notice that when 6 = the periodic solutions x = xq with wave speed c = cq of (2.6) agree 
with those of the KdV equation (2.1), which are given exphcitly by 

(2.7) xo{uj;(f>, K,k,uo) = Uo + 'i-2k^K^ cn^ (^k{uj + (j)),k), cq = uq + Sk^ k"^ — 4k^ , 

where cn(-,A;) is the Jacobi elhptic cosine function with eUiptic modulus k G [0,1) and 
K > 0,uq,(1) arc arbitrary real constants related to the Lie point symmetries of (2.1); see 
[BD]. Thus, the set of periodic traveling wave solutions of (2.1) forms a four dimensional 
manifold (3 dimensional up to translations) parameterized by uq, k, k, and (p. Note that 
such solutions are 2K{k)/K periodic, where K{k) is the complete elliptic integral of the first 
kind. 

Remark 2.2. The parameterization of the periodic traveling wave solutions of the KdV 
equation given in (2.7) is consistent with the calculations in [BD] where the authors verify the 
spectral stability of such solutions to localized perturbations using the complete integrability 
of the governing equation. However, this parameterization is not the same as that given in 
[BN], whose numerical results our analysis ultimately relies on. Indeed, in [BN] the periodic 
traveling wave solutions of (2.1) are given (up to rescaling^) as 

XsNiu; uo, q, k) = -^-^ [dn ^ , k)-^^), 

where dn(-, k) denotes the Jacobi dnoidal function with elliptic modulus k G [0, 1), and K{k) 
and E{k) denote the complete elliptic integrals of the first and second kind, respectively. 
Nevertheless, using the identity 

k^ cn2(x, k) = dn2(x, k) - {I - k^) 

we can rewrite (2.7) as 

Xo{oj) = 12k^ (dr^iniuj + <^), A;) + ^ - (1 - A;^)) , 

which, upon setting n = (]) = —ojq, and choosing uq so that 

no (i_fc2^_ m 



12^2 ^ ' K{ky 

we see that xo{oj) = Xbn{u). Thus, there is no loss of generality in choosing one parame- 
terization over the other. Furthermore, the numerical results of [BN] carry over directly to 

the cnoidal wave 'parameterization chosen here. 



''in [BN] , the authors consider the KdV equation in the form ut + &uUx + Uxxx = 0, which is equivalent 
to (2.1) via the simple rescaling u ^u. 
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Next, we consider the case < 5 ^ 1. To begin wc seek conditions guaranteeing that 
periodic traveling wave solutions of (1.1) exist for sufficiently small 6. Multiplying both 
sides by x' and rearranging, wc find that equation (2.6) may be written as 

(2.8) A + Y~'Y " + + 

hence a necessary condition for the existence of a L-periodic solution to (2.6) is 

(2.9) (xo(a;)(x'o(a;))2 - (c + l){x'o{io)f) dco = 0. 

^0 

By a straightforward computation using integration by parts and (2.6), (2.9) can be sim- 
plified to the selection principle 

(2.10) f\x'^{uj)fd^= [\x',iu)fd^, 

Jo Jo 

or, equivalently, 

pIK (k) 

/ [{cn^)']\y)dy 
Jo 

/ [{cn')"]\y)dy 
Jo 



^2K(k) 

[ 

(2-11) = =■■ F\k). 



Using the implicit function theorem, one can show that if (2.11) is satisfied, there exists a 
periodic solution of (2.6) which is 5 close to xq. As a result, we obtain a 3-dimensional 
set of periodic solutions to (1.1) parametrized by UQ,(p ^"^^ either k or n. Note that the 
limit K — >■ (i.e. A; — >■ 1) corresponds to a solitary wave and k — )• 1 (i.e. A; — >■ 0) corresponds 
to small amplitude solutions (or equivalently to the onset of the Hopf bifurcation branch). 
The above observations lead us to the following proposition. 

Proposition 2.3 ([EMR]). 5 — 0, the periodic traveling waves U{oj), u = x — ct, 
solutions of (1.1) expand (up to translations) as 



(2.12) 



J U{uj) = UQ{Ka:,uo,k,K) + 5Ui{uj) + 5'^U2{uj) + 0{5^), 
\ c = co{uQ,k,K) + 5'^C2 + 0{5^), 

where Uq^cq are defined as 
Uo{y,uo,k,K) = uo + 3k(^'^^^^ cn^ (J^y,k^ , cq = uo + {2k - 1) (^^^^^ , 

and K is determined from k via the selection principle k = Q{k) with 

f K{k)g{k) V _ 7 2{k^ -e + l)E{k) - (1 - fc^)(2 - P)i^(fc) 

V ^ ) ~ 20 (-2 + 3A:2 + SA:^ - 2k^)E{k) + [k^ + k'^- 4^ + 2)K{k) ' 
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Moreover the functions {Ui)i=i^2 o,re (respectively odd and even) solutions of the linear 
equations 

Co[Ui\ + kU'^ + K^U;;" = 0, Co[U2] +l^-C2Uoj + nU'{ + ti'Ul" = 0, 

where Cq := K^dl + dx {{Uq — Co).) is a closed linear operator acting on Lp^^{0,2K{k)) with 
densely defined domain Hp^,.{0,2K{k)). 

Proof. The explicit expansions above are determined as follows. After rescaling, continuing 
the 2K(fe)/«;-periodic wave trains of (2.1) to < (5 ^ 1 is equivalent to searching for 
2i^(A;)-periodic solutions of 

(2.13) {U - c)U' + K^U'" + 5{kU" + K^U""^ = 

for ^ > sufficiently small. We expand c, U in the limit (5 ^ as 

c = co + Sci + 0{S'^), U = UoiO + SUiiO + 0(5^). 

with Uo{uj) = xo{uj, K,k,uo) as defined in (2.7). Notice that, up to order 0(1), equation 

(2.13) is satisfied for all uo,k,K, i.e. there is no selection of a particular wave train. Now, 
identifying the 0{d) terms in (2.13) yields the equation 

(2.14) K^f/f + ((;7o - co)Ui)' - ciU'^ + kU'^ + k^U'^" = 0. 

The linear operator Co[x] = k'^x'" + {{Uq — Co)x)', defined on Hp^j.{0, 2K{k)), is Frcdholm of 
index and (1, Uq) span the kernel of its adjoint (see [BrJ, JZB] for more details). Then one 
can readily deduce that equation (2.14) has a periodic solution provided that the following 
compatibility condition is satisfied, {{Uq)'^) = K^((t^o)^) which is precisely the selection 
criterion (2.11). In order to determine ci, one has to consider higher order corrections to 
xq: in fact, ci is determined through a solvability condition on the equation for X2- This 
yields ci = (see [EMR] for more details). □ 

As a consequence, we have obtained a two dimensional manifold of (asymptotic) periodic 
solutions (identified when coinciding up to translation) parametrized by G M and wave 
number k (or alternatively the parameter k e [0, 1)). Note that the limit k — > (i.e. k 1) 
corresponds to a solitary wave and k — >■ 1 (i.e. k 0) corresponds to small amplitude 
solutions. 

3 Stability with respect to high frequency perturbations 

In this section, we begin our study of the spectral stability of periodic traveling waves 
of (1.1) in the limit S — ^ 0"*". Denote by [/(•) such a X-periodic where for notational 
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convenience, we have dropped the 6 dependence of both U and c. Linearizing (1.1) about 
U in the co-moving frame {x — ct, t) = (a;, t) leads to the hnear evolution equation 

dtv — Lv = 

governing the perturbation v of U, where L denotes the differential operator with X-periodic 
coefficients 

Lv = -{{U- c)vy - v'" - S{v" + v""). 

In the literature, there are many choices for the class of perturbations considered, each 

of which corresponds to a different domain for the above linear operator. Here, we are 
interested in perturbations of U which are spatially localized, hence we require that v{-,t) G 
L^(M) for each t > 0. Seeking separated solutions of the form v{x,t) = e^*v{x) then leads 
to the spectral ODE problem 

(3.1) Lv = \v, 
V G L2(R). 

To characterize the spectrum of the operator L, considered here as a densely defined 
operator on L^(]R), we note that as the coefficients of L are X-periodic functions of x, 
Floquet theory implies that the spectrum of L is purely continuous and that A G (t(L) if 
and only if the spectral problem (3.1) has an L°°(R) eigenfunction of the form 

(3.2) vix;X,C) = e'^Mx;X,C) 

for some ^ G [—ir/X,TT/X) and w{-) G Ll^,{0,X) Following [G, SI, S2], we find that 
substituting the ansatz (3.2) into (3.1) leads one to consider the one-parameter family of 
Bloch operators {L^}^g[_^/x,,r/x) acting on Ll^^{[0,X]) via 

(3.3) (L^w) (x) := e-^«^L [e'^wi-)] (x). 

Since the Bloch operators have compactly embedded domains Hp^^{[0,X]) in Lpgj.([0, X]), 
their spectrum consists entirely of discrete eigenvalues which, furthermore, depend contin- 
uously on the Bloch parameter ^. It follows by these standard considerations that 

(^L^R){L)= [j C^L2,,([0,X]) (^?) ; 

ee[-7r/X,7r/X) 

sec [G] for details. As a result, the spectrum of L may be decomposed into countably many 
curves such that A(^) G cr(L^) for ^ G [— tt/X, vr/X). 

The spectra A of the Bloch operators L^ may be characterized as the zero set for fixed 
^, 6 of the Evans function 

(3.4) E{X, ^, S) = det (i?(X, A, S) - e'^^Idcs^ , 
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where R{.,X) denotes the resolvent (or monodromy) matrix associated to the hnearized 
equation (3.1), that is, the solution of (L — A)i?(-,A) = and i?(0. A) = Id. Thus, the 
spectra of L consists of the union of zeros A as all values of ^ G M arc swept out. By analytic 
dependence on parameters of solutions of ODE, R, and thus E, depend analytically on all 
parameters for 5 > 0. 

In the following, we will first prove that possible unstable eigenvalues are order 0(1) + 
iO{S~^^^) by using a standard parabolic energy estimate. By a bootstrap argument based 
on an approximate diagonalisation of the first order differential system associated to (3.1), 
we show that possible unstable eigenvalues are 0(1) which implies that they are necessarily 
of order 0{S) +iO{l). We then provide an expansion in S of the Evans function as — t- in 
a bounded box close to the imaginary axis with the help of a Fenichel-type procedure and 
an iterative scheme based on the exact resolvent matrix associated to the linearized KdV 
equations. 

3.1 Boundedness of unstable eigenvalues as S ^ 

In this section, we bound the region in the unstable half plane K(A) > where the unstable 
essential spectrum of the linearized operator L may lie in the limit — > 0. Throughout, we 
use the notation \\u\\^ = \u{x)\'^dx. We begin by proving the following lemma, verifying 
that the unstable spectra is 0{S~^^^) for S sufficiently small. 

Lemma 3.1. There exist constants Ci, C2 > such that, for all S G (0, 1], the operator L 
has no L°°(M) eigenvalues with 5R(A) > Ci or (1R(A) > and ^3/^|$J(A)| > C2). 

Proof. Suppose that A is an L°° (M) eigenvalue of L and let v be a corresponding eigenfunc- 
tion. Multiplying equation (3.1) by v and integrating over one period, we obtain 



Here, we have used the fact that, by (3.2), v{x + X) = e^^-^v{x) so that l^l is X-periodic. 
Next, using the Sobolev estimate < C||?;p/2 + \\v"\\'^/{2C), valid for any C > 0, into 
the first equation yields the bound 



(3.5) 




Identifying the real and imaginary parts yields the system of equations: 



(3.6) 




(3.7) 




Letting C = 1/2 then yields 
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which verifies the stated bound on the real part of A. 

Suppose now 3ft(A) > 0. Using again the Sobolev estimate \\v'\\ < ||v||V2||^;"||i/2^ the 
imaginary part of A can be bounded as 

|5>(A)|||^f < \\U - cllooll^^f + \\v\\'/^\\v"f/\ 

Furthermore, \\v"\\ can be controlled by 

S\\v"f<{\\U'\\^ + 6)\\vf, 

which follows by setting C = 1 in (3.7) and recalling that 5R(A) > by hypothesis. Thus, 
setting K'^ = \\U'\\oo + S we deduce that 

|CJ(A)| < ||C/-c||ooifV2^-l/4+ ^3/2^-3/4^ 

which completes the proof. □ 

Remark 3.2. By slight modification, the estimates in Lemma 3.1 can be extended into the 
stable spectrum. Indeed, if C3 > then adding CsHfp to both sides of the bound (3.7) 
yields the estimate 

(5ft(A) + Cs)\\vf + 6{1 - ^)\\v"f < ^(llC/'lloo + 2C3 + 6C)\\vf. 

Thus, as long as 3f?(A) + C3 > we can repeat the proof on the estimates of imaginary parts 
to conclude 

where K3 = \\U'\\oo + 2C3 + S. 

Next, we bootstrap the estimates in Lemma 3.1 to provide a second energy estimate on 
the reduced "slow," or "KdV," block of the spectral problem (3.1). This yields a sharper 
estimate on the modulus of the possibly unstable spectrum, in particular proving that 
unstable spectra must lie in a compact region in the complex plane. Notice that this result 
relies heavily on the fact that the corresponding spectral problem for the linearized KdV 
equation about a cnoidal wave (2.7) has been explicitly solved in [BD, Sp]. 

Proposition 3.3. There exist constants Ci, C2 > such that, for all S G (0, 1], the operator 
L has no L°°(R) eigenvalues with 5R(A) > when 5R(A) > CiS or |S>(A)| > C2. 

Proof. The proof is done in two steps: first, we show that if A is an L°°(R) eigenvalue of 
L with 3?(A) > and corresponding eigenfunction v, then there exists C2 > such that 
|3'(A)| < C2. The estimate on 3fi(A) will then easily follow. To begin, let {v,\) be an 
(M)-eigenpair of (3.1) with 3f?(A) > and set x = v, y = v', z = v" + v, w = z', and 
s = c + 1, so that (3.1) may be written as the first order system 

(3.8) x' = y, y' = z — X, z' = w, dw' = —w — {U' + X)x — {U — s)y. 
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We first apply a Fenichel-type procedure and introduce wi = w + {U' + X)x + {U — s)y, 
noting then that wi satisfies 



Sw' 



{ = -wi + S(u"x + (2C/' + X)y +(U-s)iz- x)) . 

We further introduce W2 = wi — s(jJ"x + {2U' + \)y + {U — s){z — x)j so that W2 satisfies 
the equation 

6w'2 = -{l + 6\U-s))w2-6mu"' -{U-s){U' + X-6U"))x 

-^2 (^(^sU" - ([/ - sf + 6{U - s){2U' + X))y + {SU' + X + S{U- s)^){z - x)) . 

Now, by Lemma 3.1 we know that necessarily one has sft(A) + (5=^/4 |9(A) I < C for some 
constant C > 0. It follows that XS = o(l) as (5 — >■ 0, hence we may rewrite system (3.8) as 

x' = y, y' = z-x, 

(3.9) z' = W2- ([/' + A)a; - {U - s)y + s(u"x + {2U' + X)y + {U - s){z - x)) , 

5w2 = -W2 - d^x(^{z -x)-{U - s)x^ + 0{6'^{\x\ + \y\ + \z\ + |w;2|))- 

Next, we remove W2 from the equation in z by introducing the variable = z + 5w2, in 
terms of which (3.9) reads 

x' = y, y' = z^- Sw2 - X, 

4 = -iU' + X)x -{U- s)y + S(u"x + {2U' + X)y + {U - s){z^ - x)) 

(3 10) 

- (5^A((z* -x)-{U- s)x) + 0{S'^{\x\ + \y\ + \z^\ + \w2\)), 
Sw'2 = -W2 - 5'^x{{^z -x)-{U - s)x^ + 0{5'^{\x\ + \y\ + \z^\ + \w2\)). 

We further introduce the variables y = y — S'^W2, x = x and The system (3.10) 

then reads 

x' = y + 0{5'^{\x\ + \y\ + \z\ + \w2\)), y' = z + 0{5'^{\x\ + \y\ + \z\ + \w2\)), 
z' = -([/' + X)x -{U - s)y + 5(u"x + {2U' + X)v + {U - s)z\ 

(3 11) 

- 5^X{z -{U- s)x) + Oi6H\x\ + \y\ + \z\ + \w2\)), 
Sw2 = -W2 - S'^x(^z -{U- s)xj + 0{6'^{\x\ + \y\ + \z\ + \w2\)). 

In particular, by direct comparison with (2.1), we recognize the {x,y,z) equations in (3.11) 
as simply the KdV equation plus an 0(5) corrector. 

The above calculations motivate us to make a reduction to the "KdV block" of the 
spectral problem (3.1). More precisely, recalling (2.12), we write the differential system 
(3.II1, 3.II2, 3.II3) onX = {x,y,z)'^ as 

(3.12) X' = (Ao + 6{Ai + XA2) + XS^As + 0{6^)^X + 0{6^\w2\), 
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where 

(3.13) Ao 



denotes the coefficient matrix for the hnearized KdV equation, and Ai, A2, and A3 are 
defined as 

/O 0\ /000\ /OOO 

^1= ,^2 = 000 ,yl3= 00 

\ -u[ + u" -Ui + 2u' u - s J yo 1 0/ \{U-s) -1 

In order to analyze (3.12) for < 5 <C 1, we recall that in [BD] the complete integrability 
of (2.1) was used to determined a basis of solutions of X' = AqX, at least when A 7^ 0, which 
corresponds to linearized KdV equation about the periodic wave train Uq. Specifically, such 
a basis (Vi)i=i^2,3 is defined as Vi = {ui,u[,u") with Ui given by 



«.(^.A)=(l-|)exp(£ 
and r]i are solutions of the polynomial equation 



Xdy 



Uo{y)/S -c + rji 



(3.14) (r? - 46)(r? - 46)(r/ - 4^3) = A^ 

where ^1 = k"^ — 1,^2 = 2A;^ — 1,^3 = A;^. In order to deal with the limit |A| 00, we 
introduce the diagonal matrix -D(A) with 

Ai(A) = - ' ^ 



Uq/2- c + rji 



where {g{-)) denotes the average of the function g over a spatial period of U, and write a 
resolvent matrix for X' = AqX as 

R{X,co) = P{X,co)e^^^^'', 

where P{X,uj) = (Vi, V2, V3)(A, cj) is the matrix function with columns being given by the 
vector valued functions t4,i(A,a;) = e~^''''^^^^d^~^^Uk{uj),i = 1,2,3. Next we make the 
periodic change of variable X{oj) = P{X,oj)Y{oj), which is nothing but the classical change 
of variable in Floquet's theorem. In terms of Y, system (3.12) expands as 

(3.15) Y' = (l>(A) + SP-^ (^Ai + XA2 + X5A3 + 0{5)^ p)y + 0{6^ \\p-^\\ \w2\) 

as |A| — )■ 00. 

We now analyze the individual terms in (3.15) more closely. To this end, first notice 
that as |A| ^ 00 the eigenfunctions associated to the linearized KdV equation expand as 

Ui{u;,X) = (1 + (|A|-^/3\\ ^Du{x)o._ 
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It follows that as |A| ^ DO the matrix P defined above expands as 

P = f Ai A2 A3 I {1 + 0{\X\-'/^), 
\aI Ai AlJ 

where Aj := Dii{X). Thus, by a straightforward calculation we see that as |A| — >■ 00 we have 
the estimates ||P(A, •)||l°o(m) = 0(|Ap/^) and 

\\P~^{^r)\\L^{R) = 0{1), \\P~^AiP\\j^oc,^^^^)=0{l), \\P ^^3-P||L°°(R,dw) = 0(1) 

hence, using the fact that \X\S^/^ = 0(1), equation (3.15) can be rewritten as 
(3.16) Y' = (d{X) + S{XP-^A2P + 0(1)) + 0{6^^")) Y + 0{d^\w2\). 



Finally, with a near-identity change of variables of the form Y = (^Id + 0((5|A|^/^)) Y one 
can remove the non- diagonal part of (5(AP-^A2P + 0(1)) up to 0(5^/^) so that (3.16) reads 

(3.17) Y' = (l»(A) + 6diagiXP^^A2P + 0(1)) + 0(5^/^)) Y + Oi6^\w2\). 

Next, define the diagonal matrix r(A) := -D(A) + Miag(AP~-'^A2P + 0(l)) with diagonal 
entries 



r(A)ii 


= Ai + 5 1 


r^Ai(A3- 


A2) 


+ 0(|A|V3) 


r(A)22 


= A2 + 5 1 


^^A2(Ai- 


A3) 


+ 0(|A|V3) 


r(A)33 


= A3 + (5 1 


^^A3(A2- 


Ai) 


+ 0(|A|V3) 



where A := (A2 - Ai)(A3 - Ai)(A3 - A2). Prom (3.14) it follows that rji = 0(|A|2/3) 
as |A| 00, from which we see Ai(A) = 0(|A|^/^) in this limit. Introducing the polar 
coordinates A = \X\e^^'^^'^^^\ and noting that K(A) = 0(1) by Lemma 3.1, we find that 
^ = 0(|A|~^) as |A| ^00. Directly expanding the Dii{X), we have 

Ai = |A|V3e^(-/2-e/3) + o(A-V3)^ = + 0(1), A3 = fAi + 0(1), 

where j = e^'^'/3 denotes the principal third root of unity so that, in particular, we have 
the estimates 

(3.18) 5i(A,) = ^|A|V3 + o(i), sft(A3) = _^|A|V3 + o(i) 

as lAI — )■ 00. 
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With the above preparations, we are now in a position to perform the necessary energy 
estimates. Indeed, under the condition Y{x + X) = e'^'^Y{x) and W2{x + X) = e'^'^W2{x) and 
recahing that ||P(A, •)||loo(r) = 0(|Ap/3), it fohows from (3.11) that 

(3.19) ||w;2|| < C|A|<52||P(A,.)||oo(||£|| + ||y|| + \\~z\\) < CS^/'^iM + \\y\\ + \\z\\), 

where here we set Y = {x,y,z)'^. Similarly, using the bounds in (3.18), it follows from 
(3.17) that 

^5/4 

(3.20) \\y\\ + m<Cj^\\x\\. 

Inserting the bounds (3.19) and (3.20) into the x equation in (3.17) and recalling that the 
function x must be uniformly bounded on M as a function of u>, we find necessarily that 
sR(r(A)ii) = 0((5^/'') as |A| ^ oo, i.e. we have 

5R(Ai) + S (^-^ + 0(|A| V3)^ = o(<55/4) 
which, as |A| — >■ oo, reduces to 

(3.21) < -^(1 + 0{\X\-'/^)) + 6 (|A|2/3 + 0(|A|^/3)j < ^^5/4^ 

Since we have assumed 3fJ(A) > it immediately follows that |A| must indeed be bounded. 
More precisely, we deduce that there exists C2 and di > such that for all < 6 < 61, 
the operator L has no unstable eigenvalues A on L°°(R) such that |A| > C2. As we have 
already verified in Lemma 3.1 that 3fi(A) is necessarily bounded, we obtain a uniform bound 
on |S^(A)|. Moreover, it is then easy to show, by using (3.21), that, necessarily, possible 
unstable eigenvalues satisfy < 3^^(A) < CS for some constant C > 0, and the proposition 
is proved. □ 

Remark 3.4. As discussed in Remark 3.2, the estimate |$5(A)| = 0((5~3/4) ig actually valid 

so long as 5?(A) = 0(1). Thus, by repeating the argument of Proposition 3.3, one can prove 
that for any C > there exists M,6i > such that if < 6 < 5i and |A| > M, then there 
are no eigenvalues A such that 5R(A) > —C5. 

As a result of Proposition 3.3 and Remark 3.4, we have proved the following corollary. 

Corollary 3.5. Given any constant C > 0, there exist constants M,di > such that for 
< 6 < 5i we have 

(^L^R){L) C {A G C 1 3f?(A) < -CS} U {A G C I |SR(A)| < CS, \Q{X)\ < M} . 
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In summary, we have restricted the location of the unstable part of the (R)-spectrum 
of the linearized operator L to a compact subset of C, uniformly for 5 sufficiently small. 
Our next goal is to prove convergence, for a fixed ^, of the eigenvalues of the Bloch operator 

to the eigenvalues of the linearized KdV equation as (5 — > 0. This is accomplished in the 
next section through the use of the periodic Evans function. 

Remark 3.6. The structure of the argument of Proposition 3.3 may be recognized as some- 
what similar to those of arguments used in [JZ2, PZ, HLZ, BHZJ to treat other delicate 

limits in asymptotic ODE. A new aspect here is the incorporation of detailed estimates on 
the limiting system afforded by complete integrability of (KdV), which appear to be crucial 
in obtaining the final result. 

3.2 Expansion of the Evans function as 5 ^ 

In this section, we provide an expansion of both the Evans function and eigenvalues in 
the vicinity of the imaginary axis where all the eigenvalues are located at 5 = (this is 
the spectral stability result of [BD, Sp]). To this end, we will use the basis of solutions 
constructed in [BD] to build an approximation of the resolvent matrix associated to the full 
spectral problem (3.1). This leads us to the following result. 

Proposition 3.7. On any compact set X e K C C, the Evans function (3.4) of the spectral 
problem (3.1) expands, up to a nonvanishing analytic factor, as 

(3.22) E{X,^,5) = Emv(.X,0 + ^M^,0 + O{S^), < 5 « 1, 

with Ekfi^{X,^) = det (jlkdvi^^ ^) ~ e*^^/dc3^ , R^dvi-i -^) being the resolvent matrix associ- 
ated to the linearized (KdV) equation. As a consequence, for each fixed Bloch wave number 
^ G [—tt/XjTt/X) and 5 sufficiently small, if {Xs{(,))s^o £ K is an associated eigenvalue of 
then Xs{^) converges to Xq{S,), an eigenvalue of the linearized KdV equation, as 5 — > 0. 

Proof. First, we carry out a Fenichel-type computation on the spectral problem (3.1) up to 
0{5^), noting that by Corollary 3.5 the L°°(M) eigenvalues of the operator L are uniformly 
bounded in C. Recall that in the proof of Proposition 3.3 the spectral problem (3.1) was 
transformed into system (3.11): 

y + 0{6'^\w2\), y =z + 0{5'^{\x\ + \y\ + \z\ + \w2\)), 
-{U' + X)x -{U-c)y + 5 {U"x + {2U' + X)y -\- (U - s)z) 

+ 0{SWx\-h\y\ + \z\ + \w2\)), 
-W2 + 0{S'^{\x\ + \y\ + \z\ + |w;2|))- 



(3.23) 



X = 



SWn = 



Introducing Y = (x,y,z), we can thus write (3.23) as 

(3.24) Y' = A{S, X)Y -\- S'^W2 F, Sw'2 = -{1 -\- 0{6'^))w2 + S'^G'^Y, 
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where F,G e C°e^((0, X); C^) and A{5,X) = ^o(A) + 0{S) where Aq is given in (3.13). 
Further, we denote by R{-,X,5) the resolvent matrix associated to Y' = AY. It is a clear 
consequence of the regularity of the flow associated to this latter differential system that 
R{-,\,6) is regular with respect to (A, (5) and expands as R{-,X,6) = Rkdv{-,X) + 0{6) 
where Rmv is the resolvent matrix of the linearized KdV equation Y' = AqY satisfying the 
initial condition Rk(i^,{0,X) = Id^s. In order to simplify the notation in the forthcoming 
calculations, we now drop the (A, 6) dependence of resolvent matrices. 

Next, we seek to construct a basis of solutions of (3.24) valid for < 5 ^ 1. To this 
end, notice that by Duhamel's formula the system (3.24) can be equivalently written as 

Y{oj) = R{u;)Y{0) + (5M W2{'n)R{oj)R~^{r])F{ri)dri, 

Jo 

W2{u:) = exp £ tL^l^dq^ W2{0) + S £ exp ^-t^^dq^ G'^ {v)Y{r,)dv 

where here ^^'^^'^ ^ denotes a fixed analytic function of q of the specified order: for definite- 
ness, we denote this function by //(g). As a first step, we build a set of 3 eigenvectors which 
are continuations of the eigenvectors of the linearized KdV equation. For that purpose, we 
set W2{0) = and write Y as 
(3.25) 

Y{u) = R{u;)Y{Q) + 5^ £ J\^p [- j\{q)d<^ G'^ {C)Y{C)R{u)R-\^)F{n)dQdr^. 

By applying a fixed point argument in L^j,([0, X\) to (3.25), we find a set of three eigenvec- 
tors (li, W2,i)i=i,2,3 of (3.24) given by Y-i = R{^)ei + 0(6^) with aj = 6ij and W2,i = 0{5). 
To find a fourth linearly independent eigenvector of (3.24), we seek a solution (Y,W2) = 
(l4,u>2,4) such that 

(3.26) W2 = exp (^J'' ix{q)dq^ {^^^ [ ^'^'^ (" '"^^^'^^) C'^W^W^^) i 

in particular, notice then that W2{G) ^ 0. Choosing 1^(0) then so that Y{X) = gives 
(3.27) 



exp J l±^^dc^ y(a;) = -5^ j exp fi{q)dq^ R{u)R-\rj)F{r])dCdri 

+ <5^'exp p ii{q)d<^ G{OY{C)R{u)R-\rj)Fiij)dCdv. 



rX l + 0(,5^) 



We then apply a fixed point argument in weighted space e * L°°(0,X) to (3.27) to 

obtain a solution Y such that 

exp [- p ,,{q)dq^ Y{co) = -S^ £ exp ,,{q)dq^ R{OR-\'n)F{rj)dn + 0{5^). 
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Substituting this solution Y{uj) into (3.26) completes the basis of solutions of (3.24) for 
< (5 < 1. 

With the above preparations, we are now ready to expand the Evans function in S. At 
^ = 0, the resolvent matrix TZ of (3.24) reads 



7^(0,A,(5) 

whereas at ^ = 2X, it reads 



Idc3 exp (^/(f niq)dq^ 0{6^) 
exp(/,'^/x(g)dg)(l + 0(5)) 



R{X, A, S) 
0{S) 1 + 0{S) 



n{x, A, 5) = 

Therefore, it follows that 

E{X, ^,S)= det (tZ{X, A, S) - e^^^7^(0, A, 6)^ 

= -(1 + 0{S))e'^^ exp (^J^ ti{q)dq^ (det(ii;(X, A, 6) - e'^^/dcs) + 0{5^)^ , 

where we have expanded the Evans function with respect to the last column of the de- 
terminant to obtain the final equality. Recalling that R{-,X,S) = Rkdvi','^) + 0{S), the 
proposition follows. □ 

By now considering the equation E(X, ^,S) = for < S <^ 1 and applying an ap- 
propriate implicit function argument, we deduce that for each fixed ^ € [—Tr/X,ir/X) the 
eigenvalues of Bloch operator expand analytically in 5 as (5 —t- 0. 

Corollary 3.8. Let ^ G [— vr/X, tt/X) be fixed and let Xs{S,) be an eigenvalue of L^ such 
that lim^^o ^siO — ^oiC)- Then for < 5 <C 1 the eigenvalue XsiO can be expanded as 

HO = HO + SMiO + s^HO + 0(6^). 

Proof. It is a consequence of the Fcnichcl-typc reduction of (3.1) conducted above and the 
regularity of the flow of the reduced linearized problem that the Evans function is a smooth 
function of {X,^,S) on any compact subset of C x M x M.*^_. For ^ ^ 0, the eigenvalue Ao(C) 
is an isolated root of £^(-,^,0) = -Bfcd'!)(')C) so that, one has d\E{Xo{$,),^) ^ 0; see [BD] for 
more details. A straightforward application of the implicit function theorem implies that 
Xs{0 expands as Xs{0 = ^oiO + ^^iiO + ^"^^{0 + 0{S^). A similar argument holds when 
Ao(0) ^ 0. 

Now, let us consider the eigenvalue Ao(0) = of the linearized KdV equation. Notice 

then that by translation invariance and conservation of mass, coming from the conservative 
structure of (1.1), Xs{Q) = is a root of multiplicity two of E(X,Q,6) for all 6 > 0. Thus, 
the Evans function at ^ = can be expressed as E{X,0,S) = X'^{EkdvW + 0{6)) for all 
|A| sufficiently small and < 6 <^ 1, where here i?fc,j,„(A) = A~^-Bferf^(A, 0). Then A = is 
an isolated root of Ekdv with d\EkdviS^) 7^ 0, so that we can apply the implicit function 
theorem again and conclude as in the first case. □ 
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Remark 3.9. Notice that the expansion of the eigenvalues provided by Corollary 3.8 is 
precisely the one that is assumed to exist in the work of Bar and Nepomnyashchy in [BN]. 
Note, however, that this expansion is only valid for < (5 <C ^ and, moreover, is only a 
uniform asymptotic expansion for \^,\ > rj > 0, where rj is an arbitrarily small real number. 
As a result, all calculations using an expansion of the form given in Corollary 3.8 are 
valid only in this restricted, regime and, in particular, are not valid in a sufficiently small 
neighborhood of the origin in the spectral plane. 



3.3 Expansion of eigenvalues 3iS 5 ^ Q 

In the proof of Proposition 3.7 we obtained an asymptotic expansion for < J <C 1 of 
the periodic Evans function for (3.1) up to 0{5^). However, an exphcit expansion of the 

eigenvalues of such a spectral problem is often complicated to obtain by analytic Evans 
function techniques. As an alternative, here we fix a Bloch wave number ^ G [— vr/X, tt/X) 
and search directly for an expansion of the L°°(]R) eigenvalues and eigenfunctions 

u = u{-;5,^) in the form 

f m) = Ao(0 + + s'HO + o{s^) 

note that such expansions are guaranteed to exist by Corollary 3.8 and the Dunford Calculus. 
Now, recall that the spectral problem (3.1) for the operator L can be written as 



(3.29) 



/// 

u 



+ {{U - c)uy + 5{u" + u"") +Xu = 
u{x + X) = e'^^u{x), 



with u G Lpgj,([0, X]). For 5 = 0, it is known by the results of [BD] that the spectrum lies 
on the imaginary axis and it is parameterized by 

S>A = ibSVl?? - mWv - V2\\r] - %|, r/ G (-oo,?7i] U [?72,??3], 

where rji = k'^ — l,r]2 = 2k'^ — 1 and rjs = k"^ and k is the elliptic modulus associated to the 
underlying elliptic function solution ^|5_q of the KdV equation for this particular period 
X. Moreover, the Bloch wave number can be written as 



2K{k) K{k) Jo ri-k^ + dn{y,k) 



f 

Jo 



_Ntt _|_ \/|?? - mWv - V2\\'n - ml [^'^^^ dy 

r 

for some X G N. 

Before beginning our analysis of the perturbation expansion (3.28), we make some pre- 
liminary remarks concerning the spectrum of the linearized KdV operator. Let 

C:=Ll^^ = -d^,{Uo-co)-dl 

denote the linearized KdV operator, considered as a closed densely defined operator on 
L^(M), and let {C^ : ^ G [— vr/X, vr/X)} denote the associated family of Bloch operators 
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defined on Lp^^{[0, X]). By the results of [BD] , (t(£) = Wi corresponding to spectral stability 
of the underlying cnoidal wave solution Uq. Furthermore, each A G Ri \ {0} is in the 
spectrum of C with multiplicity either 1 or 3, in the sense that there exists either a unique 
^ G [— vr/X, tt/X) such that A G (t{C^) (corresponding to multiphcity 1) or else there 
exist three distinct such ^ (corresponding to multiplicity 3). Thus, when expanding such 
eigenvalues for a fixed one is essentially doing simple perturbation theory. On the other 
hand, A = is an eigenvalue of the Bloch operator jCq, corresponding to ^ = 0, with 
algebraic multiplicity three and geometric multiplicity two. Indeed, one can easily verify 
that Ker(>Co) is two dimensional and 1 G Ker(>Co); see [Br J, BrJK]. Thus, a separate 
analysis will be necessary when considering the bifurcation of the neutral modes of jCq for 
< 1. 

We now begin our perturbation analysis by considering the continuation of a fixed 
Ao G \ {0}. By above, there exists either one or three distinct Bloch wave numbers ^ 
such that Ao G a{C^). Let m(Ao) G {1,3} denote the multiplicity of Aq, as defined above, 
and let A^(Ao) = {i(,j,uoj);j = 1, . . . ,m(Ao)} denote the set of distinct Bloch wave numbers 

associated to Aq together with a corresponding function uqj G Lpgj,([0,X]) in the null- 
space of the operator jO^^ — Aq. We fix such a pair (^j,tto,j) € Ar(Ao) and set Ao(^j) = Aq, 
'^o{-',Cj) = uo,j, and insert the expansions (2.12) and (3.28) into (3.29). Collecting the 0{S'^) 
terms we find that uq must satisfy 

f + ((C/o - co)uo)' + Aono = 

which clearly holds by our choice of {Xo,uo). Continuing the expansion, identifying the 
0{S^) terms implies that ui{-;^j) and Ai(^j) must satisfy 



(3.30) 



<' + {{Uo - co)ni)' + Ao^i + Xiuo + (Uiuo)' + Uq + u'q = 
ui{x + X) = e'^^^ui{x), 



where here the function Ui is defined as in Proposition 2.3. To analyze the solvability of 
(3.30) we consider the operator /3o,iM = u'" + ((^^o ~ co)m) defined for all u G H^{Q,X) 
such that u{x + X) = e'^^^-^u{x), and note then that the operator jOqj + Aq is Fredholm 
of index on H63{0,X). In particular, we have Range(>Co,j + Aq) = Ker ((>Co,j + Aq)*)"*", 
where the adjoint operator of jCqj + Aq is given by 

{Coj + Ao)* = -dl - {Uo - co)d^ - Ao, 

defined here for all u G H^{0,X) such that u{x + X) = e^^^^u{x). Notice that if Aq / we 
easily obtain Ker(>Co,j + Ao)* via the following construction: assuming that uqj G KeT{Coj + 
Ao), one can easily verify that the function vqj^x) = J^^'^ uoj{s)ds is nontrivial, lies in 
Ker(/^o,j + Aq)*, and satisfies the boundary condition voj{x + X) = e^^^-^voj{x). Thus, for 
any such Aq and associated Bloch wave number G [— it/ X,Tr/X), we obtain a complete 
basis of Keic{Coj + Aq)*. In the case Aq = 0, corresponding to = 0, however, this 
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construction yields only constant functions: indeed, one readily finds as a consequence of 
the conservative structure of the KdV equation (2.1) that 1 G KctCq. However, we also 
note by translational invariance of (2.1) that Uq G KeiCg. In this case, we have again found 
a complete basis of KbijOq since Ker^Cp is two dimensional. 

With the above preparations, we can now obtain an explicit formula for the 0{6) correc- 
tion of the eigenvalue given in (3.28) . First, we consider the case where Aq 7^ and fix a 
such that Ao G a{C^,). Then fixing 77,0.7 ^ Kcr(/2o.j+Ao) and setting 'Uoj(.t) = f^~^'^ U()_j(s)ds 
as above, it follows by the Fredholm alternative that equation (3.30) has a solution provided 
the compatibility condition 

(3.31) (Ai-uo.i + (Uiuoj)' + iiQ j + UQ j; vqj) = 

is satisfied, where here (•, •) denotes the standard (sesquilinear) inner product on Lpgj.([0, X]). 
We now give an expression for Ai with respect to functions vqj. To this end, note by defi- 
nition we have the identity 

v'ojix) = uo,j{x + X)- uo,j{x) = {e'^^^ - l)uo,j{x), 

from which it follows that 

= ^^1^X31 X ^^.^-^^'^-^^ = 2sin(^,X/2)^ U ^^'^■^"'^■'^^J • 
Similar computations yield the following identities: 



Uiuoj)'; vo,j) = 2ze-'^^-^/2 sm{CjX/2) J Ui\uojfdx, 
*o,,) = -2e-^«^-^/2 sin(e,X/2)9 (^J u'oj^ojdx^ , 
{u'^'; vo,j) = 2e-^«^-^/2 sin(C,X/2)9 ^ u'^Ajdx^ ■ 



Taking real and imaginary parts of (3.31), assuming VQjVQjdx 7^ we can identify the 
real part and imaginary part of Ai via the relations 



(3.32) 



(/ '"kj^ojdx^ '^{Xi) = - Uilv'oj'^dx. 



Note that the 0{S) correction Ui of the underlying periodic profile U only contributes, 
up to 0(5), to the imaginary part of A. Furthermore, this contribution clearly vanishes by 
parity. Indeed, note that Iv'^J'^ = sin2(^jX/2)(l-|- |Ao ^C/qP) is an even function whereas, by 
Proposition 2.3, Ui is an odd function. As these functions are both 2X-periodic, assuming 
again that VQjVojdx ^ the integral which defines 5J(Ai) then vanishes, implying that 
^(Ai) = 0: note that this is coherent with the computations in [BN]. As a result, we have 
obtained an expansion valid up to order 0{S^) for any eigenvalue A such that A 1^=0 7^ and 
v'ojvojdx / 0. 
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Remark 3.10. In [BN], the authors numerically evaluate the expressions in (3.32) and find, 

in particular, that for any fixed Aq G Mi \ {0} that the condition Vq jVojdx ^ holds 
for each j = 1,. . . ,m(Ao). Thus, there is no loss of generality in making this assumption 
above. 

When Ao = 0, the construction of the expansion is shghtly different. Let us first remark 
that Ao = is an eigenvalue of the hnearized KdV equation which is triply covered but 
associated to the unique Floquet coefficient C = 0- Indeed the kernel of the Bloch operator 
Co is two dimensional, spanned by the functions vi = Uq and 

V2 = OmUq - ^^^dkUo = 1 - [dkCol'^dkUo . 

OkCo 

Furthermore, it is readily checked that Co{l) = —vi, hence Aq = is an eigenvalue of Co 
with algebraic multiplicity three and geometric multiplicity two; see [BrJ, BrJK] for more 
details. Similarly, we remark that when 6 ^ we have U' € KerLg, due to the translation 
invariance of (1.1), and that Lo{l) = —U'. As a result, we expect zero to be an eigenvalue 
of Lq of algebraic multiplicity two for all 6^0. Our goal now is to determine an asymptotic 
expansion of the third neutral eigenvalue of the operator Lq | = Cq for < 6 <^ 1. 

To this end, we continue the vector space spanned by {vj}j=i^2, defined above, for 
< 5 <C 1. We begin by recalling Corollary 3.8 and expanding the corresponding eigenvalues 
and eigenvectors of Lq as 

(3.33) X = dXi + 6^X2 + 0{6^), u = uq + 6ui + 6^U2 + 0{6^), 

where here we expect generically Ai 7^ 0. Substituting these expansions into the spectral 
problem Lqu = Xu, considered here on Lpgj.([0, X]), we find that collecting the 0{6^) terms 
yields -Co['"o] = 0. Thus, for some constants A^, A2 & C we can write Uq = A^vi + 742^2. 
Similarly, identifying the 0{d) terms yields the equation 

(3.34) Co[ui] + Xiuo + {Uiuo}' + Uq + Uq" = 0. 

Recalling that Cq is Fredholm of index on Hp^,^.{[0, X]) with Ker^Q = span{l,[/o}, it 
follows that equation (3.34) has a solution provided the solvability conditions 

(Al-Uo + {UllloY + Uq + Uq", 1) = 

(Aluo + {UiUoY + u'o' + u'o', Uo) = 0, 

hold. More explicitly, using the parity of uq and Uq the above solvability conditions reduce 
to 

Ai^o = 0, {Xi{v2, Uq) + {{UiV2y + v'i + Uq)) A\ = 0. 

To avoid Ai = 0, we find ^2 = ^'^d ui = Al{U[ — Ai) + A\vi + A2V2 for some constants 
yl}, ^2 G C where now we require ^4^ 7^ 0. 
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Next, we consider 0((5^) terms in (3.33) which, upon substitution into the spectral 
problem Lqu = Xu, must satisfy the equation 

(3.35) jC-o[u2] + Ai-ui + X2U0 + {Uiui + {U2 — 02)^0)' + n/ + it'l" = 

where Ui and U2 represent, respectively, the 0{S) and 0((5^) corrections of the underlying 
wave profile U: see Proposition 2.3. Using the representations of uq and ui determined 
above, equation (3.35) can be written as 

jC-o[u2] + ^? i-xj + X2U'^ + {UiU[ + {U2 - C2)u'^)' + u'C + ur) 

+ A\ {XxV2 + {UxV2)' + v'i + + A\ (X^v^ + {U^v^)' + v'{ + <') = 
or, more compactly, as 

(3.36) IIq\u2 + A\{X2 + U'^) + A\{X^ + [/()] - A?A? + {XxV2 + {JJXV2I + v'!^ + <')^2 = 0. 

Using the Fredholm alternative again, we find that equation (3.36) has a solution provided 
the solvability conditions 

-X\A\ + AiA^ = 0, (Ai(i;2, C/o) + {{UXV2I + i;^' + U^)) A\ = 

are satisfied. In particular, notice that these solvability conditions provide no requirement 
for the constant A\. Simplifying, we have thus obtained the following dispersion relation 

(3.37) A? (Ai(^;2, U^) + {{U1V2)' + v'i + v'^'- Uq)) = 

defining the 0{5) corrector Ai in (3.33). As a result Ai = is a solution and it is of 
multiplicity 2 corresponding to the Jordan block of height 2 (up to order 0(5^)). In this case, 
one has necessarily ^2 = ^ corresponding eigenfunction expands as n = A^U' + 0{6'^). 
Provided (^2, Uq) / 0, the third solution of the dispersion relation (3.37) is given by 

(3.38) ,.=_(£«):±4±^eM. 

{V2, Uo) 

In this case, A^ = XiA^ and an associated eigenvector expands as u = Uq + S(^U[ — Aq + 

X0V2) +0{6^). 

As a consequence of the above analysis, which is a rigorous version of the formal anal- 
ysis provided in [BN], for a fixed ^ G [— tt/X, tt/X) we have explicit expressions, given by 
(3.32) and (3.38), for the eigenvalues of the Bloch operator as they bifurcate from the 
eigenvalues of the associated Bloch operator for the KdV equation. In [BN], the authors 
numerically evaluate these expressions for each fixed ^ using standard elliptic function cal- 
culations: the details of these calculations are provided in Appendix A. In particular, the 
authors of [BN] find that for each fixed ^ the 0{d) correctors Ai are strictly negative for 
wave trains of (1.1) having periods lying in the interval [8.49,26.17], indicative of spectral 
stability of the associated wave trains. However, as described in Remark 3.9 this analysis is 
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only valid for < S <ti ^ and, furthermore, the expansions assumed above are only uniform 
for ^ bounded away from zero. As a result, the previous analysis is not sufficient to conclude 
spectral stability of a given periodic traveling wave U of (1.1) for any fixed 6 > 0. To make 
such a conclusion, delicate analysis in a neighborhood of the origin in the spectral plane is 
needed: this is the objective of the next section. 

Finally, we conclude this section by connecting the above analysis with the nonlinear 
stability theory developed in [BJNRZl]. As described in the introduction, nonlinear sta- 
bility of the underlying profile U under the nonlinear flow induced by (1.1) follows by 
the structural and spectral hypotheses (H1)-(H2) and (D1)-(D3). Yet as a consequence of 
the above analysis, assumptions (HI) and (D3) immediately follow from verifying that the 
0{S) corrector Ai given in (3.38) is non-zero. This observation is recorded in the following 
corollary. 

Corollary 3.11. Assume that the number Ai defined in (3.38) is non-zero. Then A = 
is a non-semi simple eigenvalue of the Bloch operator Lq: it is of algebraic multiplicity 
two and geometric multiplicity one with KerL^=o = spanjC/'} and 1 G KcrL|^Q \ KerL^=o- 
Furthermore, the return map H : M.^ U.^, {X,b,c,q) (u,u' ,u"){X,b,c,q) — b (where 
{u,u' ,u"){.,b,c,q) is solution of 

S{u' -\- u'") + u" + — -cu = q, it(0, 6, c, q) = b) 
is full rank at {X, b, c, q). 

Proof. The non-semi simplicity of the zero-eigenspace of Lq follows by the above consider- 
ations. The fact that the return map is full rank is a consequence of the fact that is of 
algebraic multiplicity 2: see [NR2, JNRZl] for more details. □ 

4 Spectrum at the origin and modulation equations 

As described above, the computations carried out in Section 3.3, which justify the formal 
approach in [BN], are only valid for |^| > > and (5 — )• 0, where > is arbitrary. 
In particular, it can be used to provide a first estimate of stability boundaries as any 
instabilities undetected would correspond to long-wavelength perturbations. Indeed it is 
verified numerically in [BN] , based on an expansion of eigenvalues similar to the one carried 
out in the previous section, that cr(L^) C {A £ Mi : Re(A) < 0} holds for all \^\ > r] > 0, rj > 
sufficiently small, for X G [Li,L2] with Li ^ 8.49 and L2 ~ 26.17: see [BN] or Appendix 
A for more details. However, one can not conclude directly to spectral stability since the 
above analysis does not rule out the presence of unstable spectrum in a sufficiently small 
neighborhood of the origin. Nevertheless, we point out that, somewhat surprisingly, these 
bounds found in [BN] are approximately those found through a direct numerical analysis 
of the spectral problem conducted recently in [BJNRZl]. In this section, we complete 
the stability analysis initiated in the previous section by studying stability of low Bloch 
numbers |^| < r] for 77 sufficiently small. In the process of verifying the spectral stability 
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hypothesis (Dl), required for the nonhnear stabihty result of [BJNRZl] to apply, we also 
prove rigorously the hypotheses (H2) ( "hyper bolicity" ) and (D2) ("dissipativity"). As a 
result, in conjunction with Corollary 3.11 and the numerical results of [BN], our results 
indicate^ the existence of nonlinear ly stable periodic traveling wave solutions of (1.1) in the 
sense of that defined in [BJNRZl]. 

s 



5 = ^1^1 Corollary 3.5 




1^1 



Figure 1: Domains in the plane where it is proved that there are no unstable eigen- 

values. The sector is described by C~^|.^| < 5 < C|^| with C ^ 1. The zone 1 corresponds 
to the domain of validity of [BN]. The zone 2 corresponds to Lemma 4.1: it is an extension 
of the analysis of [BN] and provides subcharacteristic conditions (SI), (S2), (S3) of stabil- 
ity. In zone 3, Lemma 4.3 shows that no eigenvalue can cross the imaginary axis, thus no 
unstable eigenvalue appears there. In zone 4, Lemma 4.4 proves that there are no unstable 
eigenvalues if subcharacteristic conditions (SI), (S2), (S3) are satisfied. One concludes to 
spectral stability for 6 < 6o for Sq sufficiently small. 



^Our results do not prove the existence of such solutions since it still relies on the numerical results of 
[BN]. As previously indicated, making these numerics rigorous via numerical proof would be an interesting 
direction for future investigation. 
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4.1 Spectral analysis through Evans function computations 

We begin our study of the spectrum of the hnearized operator L in a neighborhood of the 
origin by analyzing the periodic Evans function E{X,^,5) for | (A, (y)|cxRxR ^ 1- Recall 
from Proposition 3.7 that, after a suitable renormalization, the Evans function expands as 

EiX,^,S) = Ekdvi^O + SE,iX,0 + O {S\e + A')) 

for sufHciently small S > 0. Note that due to the algebraic multiplicity of the root A = of 
^kdvi'j^) = 0, the principal part of E^dv in its Taylor expansion with respect to (A, ^) is an 
homogeneous polynomial of degree 3, whereas Ei has a homogeneous polynomial of degree 
2 as a principal part in its Taylor expansion about (0,0). 

Restricting to|(A,.^)|cxR^l! the Weierstrass Preparation Theorem yields an expansion 
of the form 

EiX,^,6) = r(A,0(A - iaiiOm " i«2(0e)(A - ^3(66 + 5E,{X,0 + O {5\X^ + f)) , 

where T is an analytic function such that r(0, 0) ^ and the numbers icej{^)^ are the 
roots of the associated Evans function Ekdv(,',0 the linearized KdV equation. Using 
primarily the above asymptotic expansion of the periodic Evans function, the description 
of the spectrum of L near the origin is done in three steps. First, we show that the 
computations carried out in Section 3.3 remain valid in a sector of the form < 5 < eo|^| 
and 1^1 < rjo, for some < €0,770 ^ 1- Next, we show that eigenvalues of L^, considered 
here as a family of operators on L^^^{[0,X]) indexed by ^, can not cross the imaginary 
axis, except at ^ = 0, in a sector of the form C~^|^| < 6 < C|^| and |^| < rjc, where here 
rjc > IS small and C > 1 is arbitrary. Finally, we consider a sector |^| < eoc^ and show 
that 3?(A(^, (5)) < —6(6)^,'^ if some sub characteristic conditions are met. 

For (A,^) sufficiently small, we divide the Evans function by r(A,^) 7^ and expand it 
with respect to X,^,S as 

3 2 
(4.1) E{X, ^,6) = ll{X- ia^m + 7<5 n(A - iPU) + O {S\X' + e) + S{X' + f)) ■ 

j=l k=l 

where 7, a, 6 G M are constants, the /3k are real or complex conjugate constants, and aj(^) G 
M. Notice that since the spectral curves for the linearized KdV equation obey the symmetry 
A(— ^) = A(^), it follows that the aj are even functions of ^. Furthermore, letting ^ ^ in 
aj(^) one obtains aj(0) = o;^, J = 1, 2, 3, where the are the eigenvalues of the Whitham 
modulation system for Kortcwcg-dc Vries equation: see [BrJ, BrJK, JZB, JZl] for details. 
Concerning the a^, we note that there are well-founded numerical studies demonstrating 
that these eigenvalues are distinct for all the KdV cnoidal wave trains, i.e. that the Whitham 
modulation system for the KdV equation is strictly hyperbolic at all such solutions; see 
Section 5.1 of [BrJK] for instance. In the present case though^, we find it more appropriate 



®Note that, since all quantities are analytic, failure of strict hyperbolicity can occur in any case, if it does, 
only at isolated periods. 
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to note that the numerical results of Figure 2 in Section 4.3 below clearly shows strict 
hyperbolicity for the limiting class of KdV cnoidal waves considered here, i.e. those wave 
trains of (2.1) approachable as solutions of the KdV-KS equation (1.1) as 6 described in 
Proposition 2.3. Throughout this section, we assume that the aj are distinct and, without 
loss of generality, obey the ordering 



As a first step, we prove that the expansions carried out in the previous section, that are 
valid for |^| > rj and |A| > 77 for any rj > and 6^0 extend to a small sector < S < e\^\ 
and 1^1 < for r/, e > sufficiently small. This is the content of the following lemma. 

Lemma 4.1. There exist constants eo,?7o > and Mq > so that for allO < S < eo|^| and 

ICI < Vo, there are only three roots {Afe(^, (5)}fc=i.2,3 of E{X,(,,S) with |A| < Mq. Moreover, 
these roots are smooth functions of ^ and 5/^ and expand as 



for \ 5/(^)1 <C 1. One has A^ < 0,k = 1,2,3 if and only if the following conditions are 



(51) /3?,/30GlRand/3iV/30; 

(52) a? < < < ^5 < (once we have fixed < P^). 

(53) 7 > 0; 

Then, up to a restriction on rjo,eo, ^{Xk{^,5)) < for all < S < eo|^| and \^\ < rjQ. 

Remark 4.2. In what follows, the conditions (5*1), (52), (53) will be referred to as "the 
subcharacteristic conditions": this terminology will be justified in Section 4-2 below. There, 
we will see that the /3° are the characteristics of the first order averaged Whitham modulation 
equations for (1.1). Hence, condition (SI) above simply states that the Whitham modulation 
system for (1.1), derived for fixed 5 > Q, about the underlying wave is strictly hyperbolic. 
Note that hyperbolicity of this system, corresponding to the requirement that € M, is 
a well-known necessary condition for spectral stability to weak large-scale perturbations; 
see [Se, NR2]. We note furthermore that the condition (51) is equivalent to the spectral 
assumption {H2) necessary to invoke the nonlinear stability theory of [BJNRZl]. 

Proof. As described above, for (|A|,|^|,5) sufficiently small the equation E{\,^,5) = 
expands as 



(4.2) 




satisfied: 



3 2 



(4.3) - + - ^f^kO + o {6\\^ + e) + s{>? + e)) = 0, 



j=l k=l 
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where O ((5^(A^ + + 5(A'^ + ,^^)) stands for an analytic function of A, 5,^ of the given 
order. Now, dividing (4.3) by and setting S = 6^, X = X$, yields the equation 

3 2 

(4.4) HCX - iajiO) + 7^ - 'f^k) = O {i5\l + X') + m + ^')) • 

j=i k=i 

By comparing polynomial growth in A, it follows that there exists constants r/i, M > 0, such 
that if 1^1 + 1^1 <_r)i then |A| < M < oo. 

Now, letting ^ ^ in (4.4) one finds that necessarily 

3 

l[{X-iaj{0) = 0. 

i=i 

Since a? < < 03, the continuity of the aj{-) implies the existence of an 772 > such that 
ai{^) < 0:2(0 < o^siO for all |^| < 772- Thus, applying the implicit function theorem to 

(4.4) in a neighborhood of each aj{^) it follows that there exists three roots {Aj(^, 5)}j=i,2,3 
of (4.4), defined for ^,5 sufficiently small, which are smooth functions of ^ and 6 and can 
be expanded as 

(4.5) xM.-s) = - ,,- Ma-f;^«)-f) + om 

Ih^jWiO - «fc(C)) 

notice here we have used the fact that the aj are even functions of ^. Returning to the 
original variables via Xj{^,S) = ^Xj{^,5), we obtain the desired regularity and expansions 
for the critical eigenvalues {Xj}. 

Next, we compute the real parts of the critical eigenvalues A^ (^, 6) = ^Xj{6, (^). Recalling 
that the constants /3i,/?2 are either real or complex conjugate, as well as the fact that the 
functions ^{Xj{^,S)),aj{$,) are even function of ^, it follows by (4.5) that 

(aO-/??)(aO-/320) 

(4.6) K(A, (C, S)) = -T'^ fr (l'o7 + OiSe). 

Hence, for |^| and 6/\^\ sufficiently small, the sign of R(Xj{^,6)),j = 1,2,3 is determined 
by the sign of the real number Aj defined as 

(aO-/30)(aO-/30) 



One now needs to verify that Aj < for j = 1, 2, 3. 

For the moment, let us assume that Aj < for j = 1,2,3 and demonstrate that this 
implies the conditions (SI), (S2), and (S3) are satisfied. First, wc suppose that (3^,132 
are complex conjugates. In this case our assumption on the Aj implies that 7 7^ and 
(aO - /3?)(aO - /3^) = [a'j - /3?|2 > for each j, hence 



sgn (Aj) = -sgn (7)sgn (^("j " "fe))- 
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Since a? < 03 < Ofg, it follows that, contrary to our hypothesis, the Aj can not have all the 
same sign. Thus, it must be the case that the /3j are real and distinct, verifying condition 
(SI). Taking without loss of generality < it is now an easy computation to show that 
the signs of (^j)j=i,2,3 are the same if and only if the condition (S2) is satisfied. In this 
case, one has sgn (Aj) = — sgn (7) for j = 1, 2, 3 hence, in order to have Aj < for each j, 
it follows that condition (S3) must hold. This verifies that conditions (SI), (S2), and (S3) 
hold provided Aj < for each j = 1,2,3. Conversely, it is a elementary computation to 
show that conditions (SI), (S2), and (S3) imply Aj < for j = 1,2,3. This completes the 
proof of the lemma. □ 

In [BN], the authors computed numerically a stability index that we denote Ind(X) here. 
It is defined as follows: assume that Xn{C,S),n G N is one of the solutions of E{\,^,S) = 
expanding as A„(^, 5) = Xn{C) + ^^niO + 0(5'^), with A^j(^) an eigenvalue of the linearized 
KdV equation about a X-periodic traveling wave. Then Ind(X) is defined as 

Ind(X) = max (KA^(6)), 

?e[-7r/X,7r/X),neN 

and the authors in [BN] conclude spectral stability if Ind(X) < and instability if Ind(X) > 
0. However, the previous lemma shows that for (A, ^) sufficiently small the expansion 
assumed in the definition of A„(^, S) is apriori only valid in a sector of the form < 6 < eo|^| 
and S < 770 1^1, hence one can not conclude spectral stability from such an expansion, as is 
done in [BN], when Ind(X) < 0. However, it is easy to see that the {Aj} from Lemma 4.1 
satisfy Aj < Ind(X) for each j = 1,2,3. Hence, for the "near-KdV" wave trains of (1.1) 
satisfying Ind{X) < 0, one deduces that the subcharacteristic conditions (S1),(S2), and 
(S3) arc satisfied so that, in particular, the spectral stability in Zone 1 in Figure 4 extends 
to Zone 2 for such waves. In what follows, we prove the requirement that Ind(X) < is 
indeed sufficient to conclude to spectral stability of such a periodic wave. 

To this end, we must now deal with the range of parameters 5 > eo|^|. In what follows, 
we will assume that the conditions (S1)-(S3) are satisfied. This last part is done in two 
steps. First, we demonstrate that under this assumption no eigenvalues in a sector of the 
form C~^^ < 6 < ior C ^ 1 sufficiently large and (|A|, |^|) sufficiently small can cross 
the imaginary axis except at the degenerate boundary point {^,5) = (0,0). Then we show 
that in a sector of the form |^] < eS and {e,S) sufficiently small, one again has spectral 
stability in a neighborhood of the origin provided the subcharacteristic conditions (SI)- (S3) 
are satisfied. 

Lemma 4.3. Assume the conditions (SI), (S2), and (S3) of Lemma 4-1 hold. Then for 
all C > 1, there exists a constant rj > such that the roots of the function A 1-^ E{X,^,6) 
with \X\ < rj can not cross the imaginary axis for any \^\ < rj and ^ < (5 < C|^|, except at 
A = e = 0. 



Proof. Assume that a root A = A(^,5) of the equation E{X,^,5) = 0, with |A| sufficiently 
small so that the expansion (4.1) is valid, crosses the imaginary axis as the parameters ^ 
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and S are varied in the region 1^1 < 7/ and § < (5 < C|^|. Then there exists r = t{S,^) G M 
such that E{iT, ^,5) =0 from which, by identifying the real and imaginary part of this 
equation and using (4.1), one finds 

- "^•(^)^) = o {Sir' + e) + s\r' + e)) , 

Dividing the first of these equations by and the second one by and setting r = f ^ 
and using the estimate |^| < CS, one finds 

3 2 

- a,(0) = 0{S{f' + 1)), - P^) = 0{S{f' + 1)). 

k=i j=i 

By comparing pofynomial growth in f in the first equation it follows that |r| < M for some 
constant M. Furthermore, since the functions aj are even in ^ and |^| < CS, we have 



l[{f-a1) = 0{6), n(f-/3o^) = 0(<5). 



fc=i fc=i 
Taking 6 > sufficiently small, one finds a contradiction with condition (S2). □ 

As a consequence of Lemma 4.1 and Lemma 4.3, it follows that if the conditions (SI), 
(S2), and (S3) are satisfied, for any C sufficiently large, there exists an > such that if 
S < C\^\ and |^| < ry any solution A of E{X,^,5) = such that |A| < r? and (A,^ / (0,0) 
must satisfy 3^^(A) < 0. Next, we verify that, under the same conditions, for e > small 
enough, in the sector 5 < e\^\ and < |^| < e, all eigenvalues A, solutions of E{\,(,S) = 
such that |A| < e and (A,^) / (0,0), have a negative real part. 

Lemma 4.4. There exist constants eo,?7o > and Mq > such that for any < |^| < cqS 
and < 6 < rjo, there are only three roots {Aj(^, (^)}j=i,2,3 of E{X,^,6) with |A| < Mq. 
Moreover, assuming (SI) and 7 / 0, these roots are smooth functions of S and and 
their real part expand as ^{\^{^, 8)) = —^5 + o{5) and 

In particular, if conditions (SI), (S2), and (S3) are satisfied, then, up to possibly choosing 
i^OiVo) smaller than above, there exists 6 > such that for i = 1,2 

^iXii^,S))<~e, 0<|e|<eo<5, <5<r?o. 
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Remark 4.5. In particular, notice that Lemma 4-4 validates, under the hypothesis that 
(S1)-(S3) hold, the "dissipativity" condition (D2) of the nonlinear stability theory of [BJNRZl]. 
Furthermore, it follows that the 0{S) corrector Ai in (3.38) is non-zero provided that 7 7^ 0. 
In particular, by Corollary 3.11 it follows that the conditions (HI) and (D3) hold provided 
7 7^0. 

Proof. Recall from (4.1) that for (|A|,|^|,(5) sufficiently small the equation E{\,^,5) = 
expands as 

3 2 

(4.7) n(A - ia^) + 7^ n(A - i^U) = O {S{\^ + f) + ^^(A^ + e)) • 
i=i k=i 

Dividing this equation by 6^ and setting X = X6, = ^5 yields the equation 

3 2 

(4.8) - + 75 - = o + f ) + <^(A' + ^)) , 
j=i k=i 

which can be rewritten as 

(4.9) (1 + 0(<5))A3 + (7 + 0{5 + 0)A2 = 0(f + AO + 0(5f ) 

Hence, for |^|,<^ sufficiently small, one shows by comparing polynomial growth on A that 
|A| < M for some M > 0. Now, letting ^, 5 — ^ in (4.9) we have that 

P(A + 7) = 0, 

hence A = —7 is an isolated root of (4.8) when (^, 5) = (0, 0). Applying the implicit function 
theorem to (4.8) in a neighborhood of (A, ^,5) = (—7,0,0) implies the existence of a root 
A3 (4) ^) of (4.7), smooth in ^ and S and defined for ^, S sufficiently small, which expands as 

Furthermore, for S sufficiently small we clearly have 5?(A3(,^,5)) < 0. 

Next, we deal with the double root A = of (4.9). To this end, we assume |A| < 7/2 
and note that, from (4.9), one has 

Since |A| < 7/2, it follows that 

A2 = 0(A + P) 

hence A = 0{^) or, equivalently, A = 0{^). Now, we return to (4.7) in order to determine 
the real part of the roots \i{6,^),i = 1,2 that bifurcate from the double root A = 0. To 
this end, note that dividing (4.7) by S$,^ and setting A = A^, ^ = yields the equation 

3 2 

(4.10) e - ^"°) + 7 - '^k) = o{^s + s). 

j=i k=i 
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Now, letting ^, (5 ^ in (4.10) yields 

2 

k=l 

which by (SI) and the fact that 7 7^ has two isolate roots Ai = and A2 = iP2- 
Applying the implicit function theorem to (4.10) in a neighborhood of {\,^,6) = (Aj, 0,0) 
for j = 1,2, implies the existence of two roots {^j{^,S)}j=i^2, defined for sufficiently 
small, bifurcating from the Xj which are smooth functions of S and ^ and can be expanded 
as 

,0 , |nl=i(^«o^ 

7 J'o-Pi 



Ai(C, 6) = iPl + r "A 02 + 0{5 + 5^), 



3 



,0 , m=i^p' 



a~i 



7 Po ~ "0 



Since G M by (SI), it follows that 3fi(Aj(^, 5)), j = 1,2, are even functions of ^ that can 
be expanded as 

d2 //QO _ 0\ <:2\\-^ I'/^O _ 0^ 

(4,1) nm, s)) ^ l^iW^+o«^), uim, s)) ^ l^iW^+o«^) 

for sufficiently small. Under conditions (S1),(S2),(S3), one proves K(Aj(.^,(5)) < for 
ICI < eo*^ and < 5 < % for 770,60 sufficiently small. This completes the proof of the 
lemma □ 

In summary, it follows from Lemma 4.3 and Lemma 4.4 that no further assumptions than 
the subcharacteristic conditions (SI), (S2), and (S3) are necessary for the spectral stability 
of the underlying periodic wave train to long-wavelength perturbations. Furthermore, con- 
ditions (S1)-(S3) are validated by the numerical analysis in [BN] for all "ncar-KdV" profiles 
described in Proposition 2.3 with periods X G [8.49,26.17], demonstrating that the formal 
analysis conducted in [BN] is indeed sufficient to conclude, up to machine error the spectral 
stability of a given wave train U = Us, < S <^ 1. Finally, we also note that the analysis 
of Sections 3 and 4 demonstrate, up to machine error, that the wave trains of (1.1) which 
are found to be spectrally stable by the analysis of [BN] are nonlinearly stable, in the sense 
defined in [BJNRZl]. 

Next, we justify our terminology in referring to (S1)-(S3) as the subcharacteristic con- 
ditions by considering the formal Whitham averaged system of (1.1) about a given wave 
train U in the singular limit 5 ^ . Actually, this is precisely on the basis of the following 
formal discussion that the subcharacteristic conditions (S1)-(S3) were first conjectured to 
play a major in the small-Floquet stability of "near-KdV" waves [NR2]. 

4.2 Whitham's modulation equations 

It is now a classical result that Whitham's modulation equations for periodic waves of 
conservation laws provide an accurate description of the spectral curves at the origin, i.e. 
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of the stability of a given wave train to weak large-scale perturbations Let us mention here 
the work [Sc] in the general case, [NRl] for shallow water equations, [NR2] for KdV-KS 
either for fixed (5 > or in the KdV limit and [JZB, JZl] for the generalized Korteweg- 
de Vries equation. When considering (1.1) in the singular limit 6 — >■ 0+, however, even 
the formal derivation of such a connection is more involved. In particular, we note that 
it is not sufficient to simply let (5 —t- in the modulation equations derived for (1.1) with 
S > fixed. Instead, in this singular limit the introduction of a new set of modulation 
equations is required [NR2]. In this section, we recall the derivation of the appropriate 
modulation equations in this singular limit and emphasize in which way the previous analysis 
demonstrates their connection with the spectrum at the origin of the linearized operator 
about a given wave train. In particular, the structure of the modulation equations will 
justify our terminology referring to conditions (SI), (S2), and (S3) as the "subcharacteristic" 
conditions. 

Recall that the KdV-KS equation reads 

(4.12) dtu + a^(y ) + d^u + 6{dlu + d^u) = 0. 

Reproducing [NR2], we derive the Whitham modulation equations about a given periodic 

wave train of (4.12) in the singular limit 5 — t- O"*". To this end, we introduce the slow 
coordinates {X,T) = {ex,et), e <C 1, set 5 = 6e with 5 G (0, oo), and note that in the slow 
(X, r) variables equation (4.12) reads 

7/2 

(4.13) &ru + 5x( y) + e^d^xu + S{e^dlu + e'^d^u) = 0. 
Following [Se], we search for an expansion of u, solution of (4.13), in the form 

(4.14) u{X,T) = i7W(^^^;X,r) + £i7«(^^^,X,r) + 0(e2) 

£ £ 

with U^^\y; X,T) 1-periodic in y. Notice then that the local period of oscillation of in 
the variable y is e/dx<p, where we assume the unknown phase a priori satisfies the condition 
dx'P 7^ 0- By inserting this ansatz into (4.13) and collecting 0{e~^) terms, one finds 

(4.15) ndyU^°^ + /«c/(°)52,c/(°) + k35^C/(o) = 0, 

where = dr?!> and k = dx(t>- Equation (4.15) is recognized as the traveling wave ODE 

for the KdV equation (2.1) in the variable ny with wave speed —Q/k. As such, equation 
(4.15) has a solution provided O, = —kco{uo, k, k) where now uo, k and k are considered as 
functions of the slow variables (X, T). In this case, the solutions of (4.15) can be expressed 
as 



(4.16) 



U^^\y,X,T) = Uo{y,uo,K,k) =uo + 12k\'^cn'^{Ky,k), 
co(no, K, k) = Uo + SK^k"^ — 4k^. 
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In what follows, we derive a system of "modulation equations" describing the evolution of 
{uo, K, k) as functions of the slow variables {X, T). One such equation comes from noticing 
that the compatibility condition dxn = dx^ yields the equation 

(4.17) drK + dx{KCo{uo, K,k)) = 

for the local wave number k. 

To find other modulation equations we continue the above expansion and note that 
collecting the 0(1) terms yields an equation of the form L^dvU^^^ = . . . where L^dv is 
the operator describing the linearized evolution of the KdV equation about U^^^ . Since the 
kernel of the adjoint of Lj^dV is spanned by 1 and U^^\ solvability conditions and thus 
the needed extra eqiiations will be obtained by averaging in y against 1 and U^'^^ the 0(1) 
equation. Yet this equation being of the form 

(4.18) dr?7(°) + dx^^ = dy {■■■). 
it follows, averaging it over a single period in y, that 

(4.19) dr{Uo{-,uo,K,k)) +dx (^^{■,uo,K,k)^ =0 

must be satisfied, where here (/) := Jq f{y)dy. To obtain the other solvability condition in 
an easy way, let us first remark, following the method used in [JZl] to derive modulations 
equations for the generalized Korteweg-de Vries equation, that by multiplying (4.12) by u 
we obtain an equation of the form 

(4.20) d, + 9. (y - = 6 {{d^uf - {dluf) + dl{. • • ). 
This implies that equation (4.18) multiplied by Uq yields 

(4.21) dr (f ) + 9x (f - ^) = 5 mf - (f/^O^) + dy{- ■ ■ ). 
Averaging (4.21) over a period in y then provides the balance law 

(4.22) 0^ ) + ax (f - ^) = s - {{UHf}) . 

Together, the homogenized system (4.17,4.19,4.22) forms a closed system of three conserva- 
tion laws with a source term, called the averaged Whitham modulation system, describing 
the evolution of the quantities {uq, k, k) as functions of the slow variables {X, T). 

Again repeating [NR2], let us now comment on the previous system. As a first step in 
analyzing the modulation system (4.17,4.19,4.22), notice that the steady states are given 
by points (uq, k*, k*) G such that 

((t/^)2( • ,ns,«^r)> = ((t/^')'( fc^)> , 
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where Uq is given as in (4.16), i.e. Uq corresponds to periodic traveling waves of (4.12) 
in the hmit 5 — t- 0. Indeed, by Proposition 2.3, these are simply the cnoidal wave trains 
of the KdV equation that can be continued as solutions of (1.1). Now, letting 6^0, 
corresponding to large scale perturbations with frequency /wave number of order e S> (5, 
in the homogenized system (4.17,4.19,4.22) yields the Whitham averaged system for the 
Korteweg-de Vries equation; see [W, JZl]. As stated previously, the numerical results in 
Figure 2 in Section 4.3 below demonstrate that for all KdV cnoidal wave trains considered 
here the Whitham averaged system for (2.1) is strictly hyperbolic with eigenvalues 



Furthermore, in the limit 5 — >• oo, corresponding to a relaxation limit and large scale 
perturbations with frequency /wave number e <C (5 we obtain the relaxed system 



where here k = Q{k) is given by the selection principal in Proposition 2.3. Notice that this 
system may also be obtained directly from the Whitham averaged system of conservation 
laws for the KdV-KS equation (1.1), derived in [NR2] for fixed > as 



in the limit as 5 — t- 0. It is now well established [Se, NR2] that a necessary condition for 
spectral stability of periodic traveling waves under large scale perturbations is that system 
(4.23) be hyperbolic, i.e. have only real eigenvalues. In our analysis from Section 4.1, 
however, we assume the stronger condition that the modulation system (4.23) is strictly 
hyperbolic with eigenvalues 



this corresponds precisely to condition (SI) in Lemma 4.1. It clearly follows that in con- 
sidering only the relaxed hyperbolic system (4.23), obtained by simply letting (5 — ^ in the 
Whitham modulation equations for (1.1) derived for fixed 6 > that some information is 
lost: namely, in this particular limit we obtain no information regarding conditions (S2) 
and (S3). 

To understand the roles of conditions (S2) and (S3), we must consider rather the full 
modulation system (4.17,4.19,4.22) derived in the singular limit 5 — )• 0. For the sake of 
clarity, let us write this system with the parameterization (k, M, E) with M = {U) corre- 
sponding to the spatial average of U over a period and E = {U'^/2); see [JZB, JZl] for a 
discussion on such a parameterization of periodic wave trains of the KdV equation (2.1). 



In this parameterization, the modulation system (4.17,4.19,4.22) recovers the form ob- 
tained in [NR2] 



ai{uo, K, k) < 02(^0, K, k) < a^iuQ, k, k), V(no, k, k) G M^. 



(4.23) 





(4.24) 



Pi{uQ,k) < P2{uo,k); 



(4.25) Otk - dxi^in, M, E)) =0, drM + dxE = 0, OtE + dxQ{K, M, E) = 6R{k, M, E), 
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where n{K,M,E) = -kco{k,M,E), Q = {U^ - ZiU'^f /2) and R = {{U'^f - {U'^f). In 
the context of relaxation theory it is a classical assumption to suppose that the condition 
dER{K*,M*,E*) / is satisfied, ensuring that near the equilibrium state {k*,M*,E*) the 
equation E{k, M,E) =0 defines E implicitly in terms of (k, M); in what follows, we assume 
that this condition holds. 

Under this assumption, the subcharacteristic condition {S3) can be easily interpreted. 
Indeed, linearizing the modulation system (4.25) about the steady state (k*, M*, E*) and 
restricting to spatially homogeneous, i.e. X- independent, perturbations yields the equation 

(4.26) drK = 0, drM = 0, dxE = 6 (dsR^E + d^^MR*{K, M)) , 

where R* = R{k* , M* , E*). Considered as a constant coefficient equation in the slow 
variables {X,T), the dispersion relation of (4.26) is then given by 

(4.27) X'^{X-SdER*) = 0. 

Prom this, it is clear from our spectral analysis in Section 4.1 that the condition (53) is 
equivalent to SeR* < 0. We note that this condition is a standard assumption in the 
context of relaxation theory, and is equivalent to requiring that the manifold of solutions of 
R{k, M,E) = is stable. 

Furthermore, the dispersion relation (4.27) implies that two spectral curves bifurcate 
from the origin as one allows the period of the perturbations to vary, corresponding to 
stability or instability with respect to weak long-wavelength perturbations. It is a classical 
result [W, Yo] that a necessary condition for the stability of the steady states of (4.25) to 
such large-scale perturbations is given by the subcharacteristic condition 

(4.28) al < <a*2< < 

where here the a* and f3* denote the functions and f3j, respectively, evaluated at the 
associated steady state. Notice that in our analysis from Section 4.1, however, we assume 
the stronger condition that the inequalities in (4.28) are strict, corresponding precisely to 
condition (S2). 

In summary, we have just reviewed how conditions (SI), (S2), and (S3) were introduced 
in [NR2] as the strict subcharacteristic conditions for the relaxation type Whitham modu- 
lation system (4.25), derived from (1.1) in the singular limit 5 — > 0. As a byproduct of the 
analysis, carried out in Section 4.1, of the exact role of conditions (S1)-(S3), our present 
work have thus also rigorously validated the role of the modulation system (4.25) in the 
determination of the presence of unstable spectrum near the origin for "near KdV" waves. 
Moreover, recall that the previous analysis also implies that validation of these conditions 
follows (up to machine error) from the numerical results of [BN]. 

4.3 Numerical computation of subcharacteristic conditions 

As proved in Section 4.1 the subcharacteristic conditions (S1)-(S3) are sufficient to con- 
clude for (5 > sufficiently small the absence of unstable spectrum in a sufficiently small 
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neighborhood of the origin. Moreover, we noticed that these conditions follow directly from 
the numerical calculations of [BN]. Yet, to finish this section, we provide an independent 
verification of the conditions (S1)-(S3) using the connection to the Whitham modulation 
system of [NR2] reviewed in Section 4.2. In particular, we rely on the parameterization of 
the "near-KdV" wave trains of (1.1) described in Proposition 2.3. 

It is well known that the Whitham modulation equations for the KdV equation (2.1) 
can be diagonalized by quantities referred to as Riemann invariants; see [W]. To describe 
this diagonalization and introduce the appropriate set of Riemann invariants, we first recall 
some properties concerning the parameterization of the KdV wave trains. To begin, notice 
that traveling wave solutions of (2.1) are solutions of the form u{x,t) = u{x — ct) for some 
c G M, where the profile u{-) satisfies the equation 

uu' - cv! + v!" = 0. 

Integrating once, one finds the profile u satisfies the Hamiltonian ODE 

// 

u + — cu = a, 

for some constant of integration a G M, which can then be reduced to the form of a nonlinear 
oscillator as 

(4.29) ^ = q-W{u;a,c), W{u;a,c) = — - c— - au, 

where again q denotes a constant of integration and W represents the effective potential 
energy of the Hamiltonian ODE (4.29). On open sets of the parameter space {a,q,c) G 
the cubic polynomial q—W{u; a, c) has positive discriminant so that there exist real numbers 
ui ^ U2 < Us such that 

q - W{u;a,c) = ^{u - ui){u - U2){us - u). 

By elementary phase plane analysis, it follows that for such (a, q, c) the profile ODE (4.29) 
admits non-constant periodic solutions. Moreover, by identifying powers of u we find in this 
parameterization that 

on\ UI + U2+US I, , , . UiU2Us 

(4.30) c= , a = --{U1U2 + U1U3 + U2U3), q = — - — . 

Using straightforward elliptic integral calculations, we find that the periodic solutions of 
(4.29) can be written in terms of the Jacobi cnoidal function cn(x, k) as 

(4.31) n(0 = n2 + (n3-«2)cn^|A/^^^V^C,fel, ^ = x - ct, k" = ^^^^^. 

Us — Ui 
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In particular, notice that all solutions of (2.1) are of form (4.31) up to a Galilean shift and 

27r 

spatial translation. Letting X = — denote the period of the above wave train, it follows 

K 

again by standard elliptic function considerations that k can be expressed as 



where here 

^ ^ ^ Jo Vl — x^Vl — k'^x'^ 

denotes the complete elliptic integral of the first kind. 

Furthermore, in terms of this parameterization we note that 

,E(k) /v?\ 
(«) =ui + 2(u3-ui)-^^, {^—\=c{u)+a, 

where here (•) denotes the spatial average (in ^) over a period X and 

(4.34) E{k) = f "^^Zjfi" dx 

Jo V 1 — 

denotes the complete elliptic integral of the second kind. 

With this preparation, we can introduce the Riemann invariants {uji,uj2, for the KdV 
equation (2.1), which are defined in terms of the Ui as 

Ui +U2 Ui+ U3 U2 + U3 

In terms of this parameterization, we have 



u{^) =uji + P3-U2 + 2{u2-uJi)cn — -C,k\, 1^ = ^(j^]J — — 3 

UJ1+UJ2 + OO3 2 <^2-Wl /V , ... .E{k) 

C= , k = , {u} =LOi+U2-U3+4{U3-Ui)——, 

3 UI3 — LOi K{k) 

u'^\ 1 

c{u) + a, a = (2/3i(a;2 + ^3 - ui) + {ui + u;2 - uj3){uji + a;3 - ^2)) 



2 / ^ ' ■ ' 6 

The Whitham modulation equations for the KdV equations can be diagonalized by the 
Riemann invariants cjj, in the sense that they can be written as 

dTUJi + Vi{uJi,UJ2,UJ3)dxUJi = 0, 

where the characteristic velocities Vi are given explicitly by 

Vi{uJi,UJ2,UJ3) = 
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or, alternatively, as Vi(a;i, a;2, ws) = c + {duj^ln{K)). Clearly, the characteristic velocities 
Vi{uji,uj2, OJ^) correspond to the eigenvalues of the Whitham modulation equations for (2.1) 
about the periodic traveling wave given in (4.31) associated to (wi,W2,'^3)- To describe 
these velocities more explicitly, we find it more convenient to parameterize the problem by 
the variables wi, A = uj^ — uji and fc^ = (w2 — wi)('^3 — In terms of (wi, A,A;^), 

an elementary calculation shows that the characteristic velocities Vi can be expressed as 
14(0;!, W2, wa) = c + Ci, where Ci = ^hik) and 

k-^Kjk) k^jl - k^)K{k) (1 - k^)K{k) 

^'^^^ = E{k) - K{ky ^'^^^ = {l-k-)K{k)-E{ky ^'^^^ = E{k) ' 

with K{k), E{k) as in (4.33), (4.34) denoting elliptic integrals of the first and second kind. 
In Figure 2 we plot the characteristic velocities in terms of the period X{k) of the underlying 
KdV wave train. In particular, we see that for all k G (0, 1) the characteristic velocities are 
distinct, corresponding to satisfaction of (4.2) i.e. to strict hyperbolicity of the associated 
Whitham modulation equation. 

Next, we compute the eigenvalues of the relaxed Whitham modulation system (4.23), 
which is also the limit as (5 of the Whitham modulation system associated to (1.1) 
for fixed S > [NR2]. Recall from Proposition 2.3 that we must restrict ourselves to those 
cnoidal waves of form (4.31) such that the selection principle k = Q{k) holds. In terms 
of the {(jJi, k'^ , A) parameterization of the KdV Whitham system, this modulation system 
restricted to the "near-KdV" wave trains discussed in Proposition 2.3 can be expressed as 

(4.35) 9rK + Kdxc = 0, 9r{u) + {u)dxc + dxa = 0, 

where 

(u) = /3i + (fc^ - 1 + 4|||)A, a = -I {3P! + 2^^ + l)Ap, - {k' - ifA^) , 



and, recalling Proposition 2.3, 



^ = g{k), A(fe) = - 



3 ^ Kik)g{k) y 



To compute the eigenvalues of this relaxed modulation system, using the Galilean invariance 
of (2.1) we require {u) = 0, which is equivalent to requiring loi = —{k"^ — 1 + 4:^m^)A{k). 
This reduction thus leaves the elliptic modulus k as the only parameter of the problem. It is 
then a lengthly but straightforward calculation to show that the eigenvalues Pf{k),i = 1,2 
of (4.35) are given by the roots of the polynomial equation 



(4.36) 



A{k)X'^ - B{k)X + C{k) = 0, 
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where the coefficients are given by A{k) = G'{k), 



■6(«.-) 



Notice that the roots of (4.36) correspond to the eigenvalues /3j considered earher. In Fig- 
ure 2 we have plotted the characteristic wave speeds {a*(-^(/c))}i=i,2,3 and {(3* {X{k))}j =1^2 
as functions of the period X{k) of the underlying wave train. From these numerics, it is 
clear that the subcharacteristic conditions (SI) and (S2) are satisfied for all waves with 
period X > Xc, where the critical period is Xc ~ 8. For X < Xc, condition (S2) is violated, 
corresponding to a sideband (modulational) instability of the associated wave train. This 
threshold is consistent with the one found in [BN]. Furthermore, since the low-frequency 
stability conditions (S1)-(S3) are satisfied for all periods X > X^ we see also that the 
upper stability boundary X ~ 26.17 cannot be associated with a sideband instability, again 
consistent with the observations of [BN]. 




Figure 2: Here, we plot the characteristic velocities {aj{X{k))}j^i and {l3j{X{k))}'j^i for 
the Whitham system for Korteweg-de Vries equation and the relaxed Whitham's system 
(4.35), respectively, as functions of the period X{k) of the underlying wave train. 
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Figure 3: Here, we zoom in on the largest eigenvalues a2,3{X{k)) and l32{X(k)) in Figure 
2, noting in particular that a2{X(k)) < f32{X{k)) < a3{X{k)) for all periods X{k) > X^ of 
the underling wave train. 



Finally, we check the subcharacteristic condition (S3) and consider the spatially ho- 
mogeneous perturbations (independent of the space variable). The Whitham's equations 
read 



where u is defined by (4.31). In this setting, we use k,M = (u) and A = U3 — ui as 
parameters. One thus has 



with P{k) = 1 - k'^ + 4{k^ - 2)E{k)/K{k) + l2{E{k) / K{k)f . Next, one can show that the 
source term is written in a simpler form 



The steady states of (4.37) correspond to A = A(A;). By linearizing (4.37) about a steady 
state (A:=K, M^,, A* = A(A;=|,)) and searching for solutions that grow in time exponentially, one 



(4.37) 





R{k, M,A) = S {{{u'f) - {{u"f)) = r{k, M, A) (A(A;) - A) , 

with r{k, M, A) > and A given by: 

- _ 21 2{k^ -k^ + l)E{k) - (1 - fc^)(2 - k'^)K{k) 

^ ' ~ 20 (-2 + 3A;2 + _ 2k^)E{k) + {k^ + k^ - Ak^ + 2)K{k) ' 
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finds the dispersion relation 

A2(A- A,) = 0, 

with A* satisfying 

0^ - ) - ^) ■ 

The subcharacteristic condition (S3) is satisfied if and only if A* < 0. One clearly sees that 
this condition is independent of M^. In Figure 4, we have represented A=k = A^:/r{k^:, M^, A*) 
as a function of the period X. We clearly see that the subcharacteristic condition (S3) is 
always satisfied on the range of period [2tt, X^] with > 30. In particular, (5*3) holds 
for all near-KdV wave trains with period X > Xc, corresponding to the low-frequency 
stability boundary, and X < 26.17, corresponding to the high-frequency stability boundary 
computed in [BN]. 



-5 
-6 
-7 

m 

-9 
-10 
-11 



5 10 15 20 25 30 35 

X 

Figure 4: Here, we plot X^.{X{k)) as a function of the period X{k) of the underling wave 
train. 

Acknowledgement. Thanks to Blake Barker for the aid of his numerical Evans func- 
tion computations carried out for 5 <^ 1, some but not all of which appear in [BJNRZl]. 

A Order one eigenvalues: the computations of [BN] 

In this appendix, we supplement the analytical results of Section 3.3 by describing in our 
own notation the numerical computations carried out in [BN] determining the sign of the 
real part of the 0(6) corrector Ai(^) for a fixed ^ G [— tt/X, vr/X) for nonzero eigenvalues 
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Ao for KdV. To this end, recall that for a fixed ^ G [— vr/X, tt/X) the L°^(M) eigenvalues 
XsiO be expanded for < 5 ^ 1 as in (3.28), where we recall that for Ao(C) 7^ the 
real part of the 0{d) corrector can be found from (3.32), assuming v'q ^VQ^jdx 7^ 0, as 



(A.l) 3?(Ai) = 



notice that this is precisely formula (54) on page 593, with $0 = "Woj) for the 0{5) correction 
of non-zero KdV eigenvalues Aq found in [BN]. Using Mathematica, the authors of [BN] 
then numerically evaluate the quantity max^g[_7r/x,7r/x) 3'?(Ai(^)), which clearly must be 
non-positive to conclude stability. The details of this calculation are as follows. 
First, denote 

k' K{k)K 

Following the stability analysis for the KdV equation (2.1) presented in [Sp], the authors of 
[BN] parameterize the eigenvalues and eigenfunctions Aq and vq and the Bloch wave number 
^ as 

(A.2) vo{^)= f/"+.^^^ + "\ e-^(-+--K("), Ao = -4.'(«), ^ = 2z (c(a) - ^CM) , 

where here a and denote Weierstrass's signia- and zeta-functions, respectively, z/(z) de- 
notes the Weierstrass elliptic function with periods w = ^ and iuj' where bj' = ^^\^K(kP ~' 
Notice that ^ G M only if 3?(a) = nui, n G N. In this case, the problem is parameterized by 
a and k, since k is determined by the selection criterion n = 0{k) given by formula (34) 
on page 590 in [BN]. In [BN], the authors described the computations for a = nco + i/3 
for n = 0, 1 and /3 G [0, 2w], claiming that the other cases n > 2 do not provide any new 
results. Here, the parameter k was restricted to the interval [0, 1 — 10~^], which corresponds 
to periods X = ^ lying approximately in the interval [27r, IOtt]. In order to evaluate the 
Weierstrass elliptic functions, the usual theta functions are used: 



00 

^4/ 



Q{z) = 2^(-l)"+^gJ^"+^^ /Sin((2n - l)Trz/2K{k)) , 

n=l 

00 

ei{z) = 2^(-l)"+igJ^"+^)'/^cos((2n- l)7rz/2i^(fc)) 



n=l 



with Qo = exp(— 7rK(\/l — k'^)/K{k)). Then the various Weierstrass functions are repre- 
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sented as 

^{z) = ei + A 



@i{0)@{zVX) ) 



where ei = A = ^{u}) — v{u + iuj'). 

Using the above approach, it is found in [BN] that the quantity m.a.yi^^^_^ix^-K/X) 
is strictly negative for aU KdV wave trains with periods in the interval [8.49,26.17]. In 
particular, notice that from Figure 2 the subcharacteristic conditions (S1)-(S3) hold in this 
interval, as indicated in Section 4.1. Furthermore, the left stability boundary corresponds 
to ^ ~ 0, hence to a sideband type instability; as noted in the previous section, the right 
stability boundary does not. For each /c, and thus each period, the authors determine 
approximately the value where the functions Ai(A;,^) take their maximal values \i^rn{k). 
The points \i^m{k) = determine the boundaries of the stability region. 

As mentioned throughout our analysis, it is important to note that the analysis of [BN] 
a priori explores regions where the eigenvalues expand as 

A(<5,o = Ao(e) + <^Ai(e) + o(52), 

and is thus limited only to some particular regions of the (|C|i<^) plane. In particular, we 
stress that only the unveiling of the role of subcharacteristic conditions enables us to prove 
somewhat surprisingly that, though from the analysis of [BN] it is not possible to conclude 
spectral stability, their numerical investigation is still sufficient to complete our analysis. 

Finally, we note that another way of carrying out these computations would be to use 
instead the parameterization of eigenvalues and eigenvectors presented in [BD] . In this case, 
one has 

Ao('7) = ±8z-^|?7 - r]i\\r} - r}2\\r] - 773I, 77 e] - 00, r/i] U [772, %], 
^ _ Ntt _|_ ^\/\n - r]i\\r] - r}2\\'n - m\ f^^^^l dy 



I 

Jo 



2K{k) K{k) Jo r? - P + dn(y. A;) ' 

vo{x) = / (Ao(r/) ^jexp - / ■ — dy, 



/o Uo{z)/3- co + r]^ 

with rji = k"^ — 1, ri2 = 2k'^ — 1, rjs = k"^ and Uq the cnoidal wave given by setting k = Q{k) 
as defined in Proposition 2.3. 
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